


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1976 


Prediction of sailing boat performance. 


Munger, Edmund Colby 


Massachusetts Institute of Technology 


http://hdl.handle.net/10945/17903 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
_ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist iL Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


" \ KNOX appointed — and published -- scholarly author. 
| inp 
, LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


PREDICTION OF SAILING BOAT 
PERFORMANCE 


Edmund Colby Munger 





Sontey 





NANA 


‘ KNOX 
= <1 ws FO 


. 





7 








| 
: 
' 


PREDICTION OF SAILING BOAT PERFORMANCE 
by 
EDMUND COLBY MUNGER 
// 


B.S., United States Naval Academy 
(1967) 


Submitted in Partial Fulfillment 
of the 
Requirements for the Degree of 
NAVAL ARCHITECT 
and the Degree of 
MASTER OF SCIENCE IN NAVAL ARCHITECTURE AND MARINE ENGINEERING 
at the 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


May, 1976 





PREDICTION OF SAILING BOAT PERFORMANCE 


by 
EDMUND COLBY MUNGER 


Submitted to the Department of Ocean Engineering on May 7, 1976 
in partial fulfillment of the requirements for the degrees of 
Master of Science in Naval Architecture and Marine Engineering 
and Naval Architect. 


ABSTRACT 


Mathematical models are developed for the hydrodynamic 
and aerodynamic forces which act on a sailing boat. The 
solution of the force and moment equilibrium provides the 
Sailing boat's speed and attitude for performance prediction. 
Each force model is evaluated for its capability to make 
adequate engineering predictions. Further, expressions are 
given for the dimensions and coefficients to provide a feasible 
baseline design to which parametric variations may be made. 


The system of force equations constitute a parameterized 
system model in four degrees of freedom for steady state 
sailing conditions. The mathematical model provides a design 
tool for the naval architect conducting geometric parameter 
trade-off analysis against hydrodynamic performance. 


The principal result is a computerized model which is 
capable of predicting boat speed, heel angle, leeway angle 
and rudder angle for a geometrically defined sailing boat in 
a given wind condition. Example boats are evaluated and the 
results recorded in order to demonstrate the use and versa- 
mere y Of the model. 
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NOMENCLATURE 


keel cross-sectional area. 

projected lateral area of the hull above the LWL. 
cross-sectional area of the maximum hull station, 
or cross-sectional area of mast depending on 
context. , 

platform area of the rudder. 

total area of the sails. 

projected transverse area of the hull above the LWL. 
maximum waterline beam. 

rudder cross flow drag coefficient. 
rudder section form drag coefficient. 
Sail induced drag coefficient. 

ha tetetonvdisag COoeltricient. 

Eid weer rCELon Gdh&ag Ccoetficient. 
keel friction drag coefficient. 

hull form drag coefficient. 

mudGer Lmicel1On drag coefficient. 
rudder lift coefficient. 

Sail lift coefficient. 


rudder moment coefficient about the mean quarter 
Cnord: 


rudder normal force coefficient. 
prismatic coefficient. 

parent prismatic coefficient. 
optimum prismatic coefficient. 


rudder drag coefficient. 
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Ge} = 
LCB = 
BP opt 
LWL = 


keel side force coefficient 
slender-body side force coefficient. 
rudder side force coefficient. 

sail drag 

offwind sail drag. 

aerodynamic hull force. 

Froude Number. 

maximum canoe body draft. 

heeling moment. 

hydrostatic righting moment 

keel heeling moment. 

rudder heeling moment. 

Sail heeling moment. 

load waterline length. 

Sau le lies 

CHiwind Sale Lift. 

position of center of buoyancy (%SLWL). 
parent position of center of buoyancy. 
optimum location of center of buoyancy. 
load waterline length. 

pitching moment 

yawing moment 

keel yawing moment. 

rudder yawing moment. 

slender-body yawing moment. 

yawing moment due to sail driving force. 


yawing moment due to sail side force. 
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vector of assumed variable values. 
=assumed value of individual variable. 
total frictional resistance. 

total form resistance. 

Reynold's Number. 

upright resistance. 

wave resistance. 

parent wave resistance. 

vector of force equilibrium errors. 

= individual force equilibrium errors. 
fairbody wetted surface. 

keel wetted surface. 

rudder wetted surface. 

canoe body maximum draft. 

boat velocity in the Xo direction. 
apparent wind velocity. 

boat velocity parellel to the y axis. 
true wind velocity. 

force in the XS @irection. 

hull resistance. 

keel and hull induced drag. 


axis on the plane of the free surface in the 
GiaBecerTOnwo: boat velocity. 


rudder induced drag. 
Sime araving LOrce. 
force in the — Grreection,; Side force. 


keel side force. 
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an axis on the plane of the free surface at 
right angles to boat velocity. 


rudder side force. 
Sail side force. 
force in the Zo direction, heave force. 


an ax1sS perpendicular to the free surface 
downward. 


. fortriangle aspect ratio. 


mainsail aspect ratio. 

rudder aspect ratio. 

total sail area aspect ratio. 

draft at forward perpendicular. 

maximum draft. 

rudder span. 

intermediate draft measurements. 

rudder root chord. 

rudder tip chord. 

rudder quarter chord. 

keel mean chord. 

rudder mean chord. 

metacentric height/LWL. 

non-linear hydrostatic righting moment coefficient. 
Froude Number dependent righting moment coefficient. 
spline cubic function. 

acceleration due to gravity. 

canoe body draft at maximum draft measurement. 
Sinkage. 
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t = taper ratio. 


r 
Xx = axis in the load waterline plane along the 
Gemeerline of the boat. 
xy = Heorwzontcal Gistance to the maximum draft 
measurement. 
> horizontal distance to the aerodynamic center 
ole Gin cttor-e On the hull. 
oat = hO@uiwemmedal distance to center of effort on 
P the keel. 
a = hemmAomtel diStance to the center of effort on 
P the rudder. 
ee horizontal distance to center of effort on 
aA the sail. 
Xo = main boom length. 
Be = base of the fortriangle. 
cts = horizontal distance to the mast. 
ie = U0 location of maximum draft. 
Xx es horizontal distance from the rudder mean chord 
P quarter point to the rudder center of effort. 
me = horizontal distance to the rudder stock. 
xy 2.37 horizental distance from maximum draft to 
ad intermediate draft measurements. 
y = axis orthaginal to the x axis to starboard 
in the waterline plane. 
Z = axis orthaginal to the waterline plane, downward. 
Z = vertical distance to the center of effort on 
cpk 
the keel. 
Z = vertical distance to the center of effort on 
Sore 
the rudder. 
Zcpxs = vertical distance to the sail driving force. 
Zopys = Veterledimatstance to the sail side force. 
Zep = freeboard. 


Ie 





vertical distance to the center of gravity. 
mast height. 
MacdmynOLSt.. 


vertical distance from the rudder root chord to 
the center of effort on the rudder. 


vertical distance to the rudder root chord. 
volumetric displacement of the fairbody. 
volumetric displacement of the keel. 

weight of the boat in tons. 

rudder angle. 

apparent wind angle. 

upright apparent wind angle. 

true wind angle. 

variable increment. 

pitch angle. 


angle between boat centerline and aerodynamic 
hull force. 


kinematic viscosity of water. 
Bs dA 5 0 

density of air. 

density of water. 

heel angle. 

leeway angle. 


rudder sweep angle. 


1 





fee RODUCTION 

1.1 General 

All designers of engineering systems desire a procedure 
for predicting the performance of their designs prior to 
meme CcOnStruction. Until lately, prediction of a sailing 
yacht's performance has been limited to evaluation of scale 
models and intuition. The development of the high speed 
computer and numerical techniques for rapid solution of non- 
linear equations has made the use of mathematical models a 
reasonable alternative for performance prediction. It is 
the purpose of this paper to present mathematical models for 
the forces and moments acting on a sailing boat. The solution 
of the force and moment equilibrium will provide the boat's 
speed and attitude for performance prediction. 

Fach force is parameterized in dimensions and coefficients, 
generally utilized by designers, to make the model a useful 
design tool for the naval architect. Further, expressions 
are given for the required dimensions and coefficients to 
provide a feasible baseline design from which parametric 
variations may be made. Each force model is evaluated for 
1ts capability to make adequate engineering DredilceELoOns, 


iMelieaing inherent limitations. 


Paces eomical Development 


Prom the earliest days, yacht design has been called an 
rt see thneugn the use Of wind propulsion for ships is ancient, 


the action of the forces and moments on the sailing ship is 
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complex and multidimensional. Historically there has been 
no technique for simple evaluation of the effect a change in 
one parameter has on the forces acting on the boat. 

Until fifty years ago deSign evolved slowly through ober- 
vation of existing ships. Due to the pebondh iy to evaluate 
a change prior to Meerooraticn into a newly constructed 
yacht, new designs tended to consist of only conservative 
changes to a previously successful design. This was under- 
standable for failure of an innovation in the final yacht 
could be professionally catastrophic for the designer. 

In the early 1900's, techniques were developed for pre- 
dicting yacht performance from model tests. A deSigner could 
now build a relatively inexpensive model and evaluate it 
against other model results prior to incorporating a modifica- 
tion into a new deSign. Adequate prediction of full scale 
results have been limited by adverse scale effects, but the 
usefulness of these tests is unquestioned. Unfortunately, 
the cost of model testing still precludes extenSive parametric 
trade-offs. Each parametric change requires construction and 
testing of individual models. Usually, the limited design 
budget for one yacht precludes any systematic variation of 
parameters, and testing is limited to one or two candidate 
models. The appeal of a parameterized mathematical model 


which is inexpensive to utilize iS apparent. 
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The idea of using mathematical expressions for describing 
yacht performance is by no means a new concept. The rating 
rules used to handicap yachts are an attempt to apply 
mathematical models to boat performance. Because rating 
rules attempt to measure a boat's speed potential for a large 
range of conditions and reduce the results to a single number, 
it is easily understood why they have little application to 
predicting performance for a given wind direction and speed. 

It was not until the second half of this century that 
towing tank test results were combined with sail force 
coefficients in a computer to solve force and moment equili- 
briums. These programs could provide prediction of performance 
for a given condition within the limitations of the data. 
These solutions were for the specific boat tested and a 
scaled family of geometrically similar boats. 

In the early 1970's attempts for steady state solution 
from mathematical models began to be published. One pioneer 
in the field has been John S. Letcher, publishing parts of 
his model in references (1), (2), and (3). Further work has 
been done by Myers (4) and Kerwin (5). The computational 
technique for solution of the force equilibrium condition 
has been described by Hazen (6). 

By parameterizing forces in four degrees of freedom, the 


model developed in this paper grows from these previous works. 
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1.3 Overview 

In order to resolve the aerodynamic, hydrodynamic and 
Pyamostatic forces acting on a sailing boat a coordinate 
system and origin are established. Achieving an equilibrium 
Semeereces and moments requires that the summation of all 
forces and moments acting either through the origin and about 
the origin are equal to zero. To facilitate modeling, 
component forces which act independently or nearly so are 
identified and examined. Once the forces acting on the system 
and the location and direction through which they act are 
identified, each may be modeled and superimposed to give 
the required equilibrium solution. 

1.3.1 Coordinate System 

Figure 1 presents a coordinate system which describes 
the sailing boats velocity and attitude. The coordinate 
system designated by the capital letters Kor You and Z, are 
located on the free surface with Xo and Yo axis in the 
horizontal plane. The oe axis lies in the direction of the 
boat's motion and the Yo axis is to starboard at right angles 
eOments velocity. The Zo axis is vertically downward. 

A second coordinate system is fixed in the boat at the 
forward end of the waterline as presented in Figure 2. This 
coordinate system is designated by the small letters x, y, 
and z. The x axis is oriented along the fore and aft center- 


line on the boat's waterline with the y axis in the transverse 


ey 





Xe 


AZ 
; Al 


Zo 2 


Figure | - Coordinate Systems 
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Boat Coordinate System 
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Girection to starboard. The z axis desends orthogonally to 
the x-y plane in the direction of the keel. 

The angular differences between the two coordinate 
systems describe the attitude of the boat as it moves through 
the water. The positive direction for these anales are also 
presented in Figure 2. Their direction was selected to Puo- 
vide positive values when the boat is in equilibrium under 
normal sailing conditions with a weather helm. The order in 
which these angles are taken is important to the final atti- 
tude. The boat is first rotated about the a axis through 
the angle w. The angle ~ is identical to the boat's leeway 
angle. The boat is then rotated about the X axis, an angle 9, 
which is the boat's heel angle. The boat is then rotated 
about the Yo axis, an angle 6. 8 is the angle by 
which the boat pitches down at the bow and is also the angle 
of the axis about which the boat heels. Finally the angle ort 
is the angle the rudder is deflected from the boat's centerline. 

A sailing boat responding to the forces acting upon it 
and conforming hydrostatically to deformations in the free 
surface would heave in the positive Zi direction an amount s. 

The selection of these coordinate systems allows the 
modeling of forces in either coordinate system. When modeling 
velocity related forces, such as wave drag and lift, the 
coordinate system fixed in the water is more attractive 
computationally. On the other hand, hydrostatic forces are 


more easily computed in the coordinate system fixed in the 
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boat. In the final solution of force equilibrium, forces 
and moments are resolved in the Xoe Yor and Zo system except 
the heeling moment which is resolved about the x axis. 
It 1S a common practice in naval architecture to describe 
a vessel's dimensions and coefficients in a positive sense. 
With few exceptions this is done throughout this paper. 
The positive direction of forces and moments is defined in 
Seerion 1.3.2. 
eee Lei brium Of Forces 
The forces and moments acting upon a Sailing boat 
in these coordinate systems may be divided into six components: 
X = Forces in the positive Xo direction 


Y = Forces in the positive Xo dihBecElon 


N 
il 


Forces in the positive Zo direction 

K = Moments about the x axis 

M = Moments about the y axis 

N = Moments about the Zo ax1s 
The positive direction of the moments K, M, and N are taken 
in the same sense as the angles ¢, 6, w respectively as 
indicated in Figure 2. Under normal conditions sail aero- 
dynamic forces will create positive forces and moments on 
the boat. 

Through out the remainder of this paper the moments K, 


M, and N may be referred to interchangeably as forces or 


moments. 


Zen 








The resolution of forces into equilibrium requires that 
the summation of all the forces acting in the six directions 


be equal to zero, i.e.; 


eo tt 


where "i" designates individual aerodynamic or hydrodynamic 
forces or moments acting in the prescribed direction. An 
example of a contributing force is Kor Enewscad SaGniving force. 
This is one of a number of forces acting along the Xo axis. 

The achievement of an equilibrium of forces in the six 
directions requires the solution of six non-linear simultaneous 
equations. To obtain a unique solution, six variables are 


required. A selection of variables for the six degrees of 


freedom problem is: 


(U,, eo) fae a) 


The only unidentified variable is U which is the boat's 


b! 
velocity in the positive Xo direction. Each variable contri- 
butes a major effect to a corresponding force. The boat's 
resistance is strongly dependent on boat speed, UL: As lee- 
way angle, w, strongly contributes to side force, so does 


heel angle, 6, to the righting moment and rudder angle, Qs xe) 


the yaw balance. 


Z2 
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The present model excludes the solution in heave and 
Beech. This exclusion of two degrees of freedom was not due 
to their lack of importance, but the lack of an adequate 
parametric model. In reality the axis about which the boat 
heels may be pitched as severely as 15°. This pitch makes 
a measurable change in the boat's angle of incidence when 
heeled to large angles. Further, due to the drop in dynamic 
pressure as the water passes around the hull : certain amount 
of sinkage will occur as a measure of the hydrostatic 
reaction. Besides the effect on form drag, this dynamic 
pressure and resulting sSinkage also increases the boat's 
wetted surface and viscous resistance. The magnitude of 
these effects is small enough that they may be reasonably 
ignored. In the generation of mathematical models for the 
remaining forces, emperical data was used that came from 
experimental tests which included the effects of both 
beech and heave. 

The present model fixes the pitch angle, 9, and sinkage, 


Ss, at zero. The variables are reduced to the following four: 


( Une Vr Os a) 


The remaining major forces acting on a sailing boat are 
resolved into X, Y, K, and N components and are added together 


to give the following force equilibriums: 


Za 





UX. 
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BAY 


ZK. 


BN 


, thy 56 ee 
Wert ove YY. = 0 
S r 
Kae KK 
S 


ort kK. = 0 


N ae at Zon. +N. = 0 


The individual forces are described as follows: 


Satiedrivangetorce 


Hull resistance 


Hull and keel induced drag 


Rudder induced drag 


Sail side force 


Hull and keel side force 


Rudde~ side force 


Sail heeling moment 


Hobe entang moment 


Keel heeling moment 


Rudder heeling moment 
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Ne = Yawing moment from sail driving force 
ae = Yawing moment from sail side force 
Ny = Yawing moment from keel side force 
a = Yawing moment from rudder side force. 


Each individual force is defined in later sections as a 
function of the boat's geometry, a true wind vector and the 


four variables U Ww, o, and a. The solution of equations 


b’ 
1 through 4 for a given boat geometry and wind condition 
provides values for boat speed and attitude. From these 


values the hydrodynamic performance of the specific sailing 


boat may then be evaluated. 
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II. FORCE AND MOMENT MODELS 

2.1 Hydrodynamic Forces 

The hydrodynamic forces acting on a sailing boat may be 
conveniently divided into four areas for modeling. The four 
feetomace hull forces, keel™forces, rudder forces, and hydro- 
static forces. Hull forces consist of resistance from wave 
making, viscous effects, and form drag. Keel forces are 
made up of lift on the keel and hull, and drag induced by 
this lift. Rudder forces in the same way consist of lift 
and induced drag from the rudder. Finally, hydrostatic 
forces are caused by the change in position of the center 
of buoyancy in the hull. Each of these areas is presented 
individually and models are developed for their contribution 
to the total force equilibrium. 

2.1.1 Hull Resistance 
The estimation of hull resistance is a problem faced 

by naval architects in the design of every ocean going 
vehicle. The difference, in predicting a sailing boat's 
resistance, is in the determination of the effects changes in 
attitude have on these calculations. The approach taken is 
to determine a model for up-right resistance and modify this 
resistance as necessary to account for changes in the boat's 


attitude. 
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Besiohnt Resistance 


The technique commonly used for prediction of resistance 
from tow tank data is based on Froude's Hypothesis. Basically, 
the hypothesis states that the total upright resistance may 
be divided into two components. The first is frictional or 
viscous resistance. Frictional resistance is dependent on 


Reynold's Number, a non-dimensional parameter defined as: 


R. _ UDE/ 4, 0 
2 
The second component is residual resistance. The residual 
resistance is the combination of drag due to the generation 
of waves and the effects of eddy making and form drag. The 
effect of the hull generating waves is a gravitational effect, 
largely dependent on Froude Number. This non-dimensional 


coefficient is defined as: 


Though Froude left the separation of the elements of 
resistance in these two components, form drag and eddy making 
resistance have little dependence on either Froude or Reynolds 
Numbers over the velocity ranges encountered by sailing 
boats. For this reason it is logical to further separate 


residual resistance into wave resistance and form resistance. 
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The separate components of upright resistance are assumed 
independent of each other so that the total resistance of a 
boat at a given speed may be estimated by superposition of 


the forces. Total upright resistance is then: 


Upright resistance = frictional resistance 


+ wave resistance + form resistance. 


Each component of upright resistance is examined in detail 
and mathematical expressions are developed to quantify their 
magnitude. 

Wave Resistance 

The mathematical model for wave resistance is based on 
a generalized form of Michell's Integral. This expression 
is evaluated for a base hull form and is corrected for 
parametric deviation from this parent form. This expression 
is emperically adjusted to give accurate prediction of the 
towing tank results of Antiope (3). 

Michell (7) postulated that wave resistance for thin 
ships depends upon the computation of pressure changes, O57 
due to the wave pattern. His expression for this resistance 


1s 
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where 


2 2 
ma f ff ons e A g2/U,, ae CLES ax dz 6 
Ox Zz 
-e Q U 
b 
e K(x) 2 2 
w= | f 23 A gz/U;y, 299 (SES) Gace I 
Ox 2 
-e 0 Uy, 


Another form of these equations was derived by Weiblum (8). 


A general expression of his result is 
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By assuming that the block coefficient, Che varies very little 
between yacht forms, equation 8 may be written as: 
R. = _ (a c(i, form coefficients, Fy) 9 


Wg “ (Z/100) 2 ae 


The number of experiments conducted to validate Michell's 
result is small. Wigley (9) compared model results to theory 
and observed the interference humps in the resistance curve 
were overpredicted by theory. Further comparison between 
vessels of constant draft and varying beam done by Wigley (10), 
Weinblum (11), and Mumford (12) indicate for slender ships 
wave resistance is approximately proportional to the square 
of the beam as estimated by Michell's Integral, but as the 
beam to length increases dependence of resistance on beam 


appears to be less than quadratic. It is expected that the 
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high beam to length ratios encountered in sailing boats will 
strain the slender ship limits of Michell's Integral. 

The basic form of the wave resistance iS rewritten as: 

L3 A n ae 

R = Oe | ae E(L/TO, fonm COcEEICTenEeS, EF) 1L10 
where n 1s between one and two. The actual value is 
determined by emperical data discussed in the following 
sections. < 

The identification of the most important parameters 
affecting wave resistance is required. Todd (13) dida 
regression analysis of over two hundred destroyers and found 
the following coefficients had the highest correlation between 


the speed length ratios of 0.8 and 1.25. 
3 
(77a 100)", oy Les; L/T F) 


For geometrically similar boats we know directly from 
Froude's Hypothesis that the wave resistance scales is the 
cube of the characteristic length. The general form of 
Mache lis integral, equation 10, is used for the dependence 
on displacement to length ratio. By assuming base value for 
DisiSMatic coefficient, Cao! longitudinal position of the 
center of buoyancy, LCBO: anemmengeh to draft ratio, L/T Gi 
an expression for the wave resistance of a parent form may 


be written: 
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ae = £, (U,/VvL) L° [—4——~]" dk 
Wwe OS oin/ZaLieien 


eels, LCB., and L/T 


po on 


To predict the wave resistance of a sailing boat, the 
resistance of the parent boat is adjusted by correction 
factors based on the boat's geometric deviation from the 


parent. The wave resistance equation is 


La Rk PeetrCe GOrrection factor + LCB correction 
W Wp p 


factor + L/T SGEheCeLIOn Laccton) 2 


Fach component of this equation is defined and discussed 
in the following sections. 

Characteristic Length - L 

The characteristic hydrodynamic length or "effective 
length" of a sailing boat is difficult to quantify. This 
has been a difficulty in the maintenance of all offshore rating 
rules as demonstrated by the extremes to which designers have 
gone to make rated length look shorter than the boat's actual 
hydrodynamic length. For the purposes of this model, L is 
determined to be the load water line, LWL. This creates 
little difficulty when examining the comparitive performance 


of similar forms, but subjective adjustment would be required 


Sy 





when comparing plumb bow, transom stern boats to long ended 
boats with shallow buttocks and entrance profiles. To some 
extent the latter might increase in length with Froude Number 
as the boat will be supported more by her overhangs. 


Parent Froude Number Dependence - £, (U,/VL) 


The speed length ratio, U/VvL, was adopted instead of 
the dimensionally correct variable F_ = U.//gL due tomthe 
common usage and the comfortable feel most yacht designers 
Mave with) speed length ratio. The function, £, (U,/VvL) aq justs 
the wave resistance equation to accurately predict the 
resistance of the yacht Antiope. The function was determined 
by inputting Antiope's geometry into the model and selecting 
£, (U,/vL) to insure proper correlation with the full scale 
towing test results reported in (3) and (14). Values of the 
function were eStablished at each data point given in (14) 
and a spline cubic curve was fitted to give the curve in 


Broure 3 for n = 2.0. 


n 


Displacement to Length Factor - [A/(L/100) >} 


The high beam to length ratios, B/L, encountered in 
sailing boats requires that the exponent, n, be determined 
from emperical data. Unfortunately, the data base from tests 
of actual sailing boats is small. The scale effects pointed 
out in reference (15) precludes the use of results from towing 


tests of small scale yacht forms. 
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Figure 3 for n = 2.0. 
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Becplacement to Length Factor —- [A/(L/100) 2] 


The high beam to length ratios, B/L, encountered in 
sailing boats requires that the exponent, n, be determined 
from emperical data. Unfortunately, the data base from tests 
of actual sailing boats is small. The scale effects pointed 
out in reference (15) precludes the use of results from towing 


tests of small scale yacht forms. 
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Examination of Taylor Standard Series (16) at prismatic 
coefficients typical in yacht design indicate that a value 
of n between 1.8 and 2.0 appears reasonable. Though the 
prismatic coefficients, length to draft and displacement to 
length ratios in reference (16) are typical of sailing boats, 
the block coefficient, Chr is high and not characteristic OF 
yacht forms. 

Figure 4 is a plot of two curves predicting the wave 
resistance of a Standfast 43 sailboat by the mathematical 
model with n equal one and two. Also plotted are the data 
points from the towing tests of a one sixth scale model 
conducted by M.I.T. and Delft University. Figure 5 compares 
the full scale test (15) of the one ton yacht Bullet to two 
theoretical curves. The form drag data from both yachts was 
taken as the residuary drag at a speed length ratio equal to 
0.4 and subtracted from the residuary resistance versus speed 
curve to obtain the wave drag. 

Two iS maintained as the value of n in the computerized 
model appended to this paper. Though comparisons made in 
Figures 4 and 5 indicate that n = 2 may slightly overestimate 
the effect displacement to length ratio has on wave making 
resistance, the limited sample size available for correlation 


does not justify adjustment from the theoretical value. 
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@eGorrecteion Factor 
cal al ae 


The prismatic correction factor corrects the magnitude 
of the sailing boat's wave resistance for deviations away 
from the base prismatic coefficient, ook The base prismatic 


coefficient, C30! was selected to be the optimum prismatic 


coefficient, C for each speed length ratio. Taylor's 


popt’ 
Standard Series (16) data are the only tests available of 
large models where prismatic coefficient was systematically 
varied. Examination of these data reveals a definite optimum 


Prismatic coefficient, C for each tabulated speed length 


jaKejoLe 
Beto. Figure 6 is a spline cubic fit of this optimum 
prismatic taken at the displacement to length ratio equal to 
250 and plotted versus speed length ratio. This optimum 


Beismatic coefficient, C as a function of the speed 


popt’ 
length ratio is selected as the base prismatic coefficient, 
Coo’ of the parent equation. This is done so that deviations 
from this optimum will always result in an increase in wave 
making resistance. 

From Taylor's Standard Series, reference (16), the 
percentage increase in resistance for deviation from the 
optimum prismatic coefficient was tabulated for each speed 
length ratio. From this emperical data the expression for 
the S3 correction factor is generated. 


Cy eormection factor = £. (U,/VL) |100(C “cy |e! 13 


popt 


or, 





meeure / 1S a plot of the spline cubic fit of the correction 
factor's dependence on speed length ratio, £,(U,/VL) . A 
reasonable definition for the sailing boat's prismatic 
coefficient, ou must be established to insure quantitative 
agreement between the Taylor Series, reference (16), forms and 
the yacht forms. A method for weighting the volume and 
sectional area of the keel in the calculation of the longi- 
tudinal prismatic coefficient is developed. Equations 6 and 

7 of Michell's Integral indicate that displacement should be 


weighted exponentially in the vertical direction bv 


a7\g2/U,. 

There are two common practices in calculating the oe OE 
a sailing boat. The first is to include the total volume 
and sectional area of the keel into the calculation. The 
inclusion of all the deeply placed sectional area of the keel 
would overestimate its effect on the wave making resistance 
when compared to the keel's actual effect on the wave train. 
On the other hand, the second method is to treat the keel 
as an appendage and calculate only a fairbody or canoe 
prismatic. This ignores the keel's displacement effect 
altogether. Typical values of a for the same yacht calculated 
both ways are 0.52 and 0.58 with and without the keel 
included respectively. If the Taylor Series, reference (16), 
optimum is to be used, the corresponding prismatic coefficient 


would be somewhere in between these two values. 
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One would expect the volume of the keel to have its 
largest effect when the sailing boat is sitting on a wave 
of approximately the same length as the waterline. At this 
point the wave Resusk would bring the keel the closest to 
the free surface. An estimation of keel's proper effect on 
prismatic coefficient, > may be made by weighting the keel's 
area and volume by gol ee where ee 1s the vertical Froude 


number evaluated at speed length ratio equal to one. The 


longitudinal prismatic coefficient, or is defined as: 


V+ Ves 
* 
(A, aE Ax) L 
where 
mopNa (12689) 
Ve = V, e ° ‘ 15 
Kk k 
and 
2 
me e 32-2d/ (1.6891) 16 
ins k 
U- and ey are the displaced volume of the canoe body and 
keel respectively. An and A, are the maximum cross sectional 


area of the canoe body and the cross sectional area of the 
keel at that station respectively. The depth, d, is measured 
from the static waterline to the center of the keel's displaced 


volume. 
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L/T Correction Factor 





Taylor Standard Series (16) tabulates data on ships at 
Jength to draft ratios equal to eight and sixteen. Taylor 
Series (16) asserted that linear interpolation between these 
ratios was adequate. 

The base length to draft ratio, L/T oo? Gis o. Ss 
selected from Taylor Series (16) data and the correction 


factor is 
oD correction factor = ~£,(U,/VL) (16-L/T) 7 


where £,(U,/vL) is a spline cubic fit of the speed length 
effect as shown by Taylor (16). (£,(U,//L) Hem ploreecdmsan 
Figure 8. 

The major differences between the ships tested for 
Taylor's Standard Series (16) and yacht forms are in the 
Shape of the sections. The Taylor form typically has a 
maximum section coefficient, Ca of 0.98 while sailing boats 
range between 0.60 and 0.70. As a result, the vertical 
distribution of volume between the two forms differ to a 
large extent. It is not completely accurate to assume 
that vertical shifts in the volumetric distribution of the 
two forms will have a similar effect on wave making resistance. 
The model utilizes the best data available to date, but the 
need exists for a large scale yacht series examining just 


such a systematic parameter variation. 
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mem Correction Factor 

In the same manner as the Ss correction factor, the LCB 
correction factor is based on the actual LCB and an optimum 
LCB position. A summation of recommended optimum LCB positions 
Pema function of block coefficient and prismatic coefficient 
was presented by Todd (17). There is reasonable agreement 
among the recommended curves with less than 1.5% variation 
between extreme curves in the area of interest. Comparing 
the LCB curves with common practice in yacht design a curve 
attributed to Professor Troost in (17) was adopted. In 


meaitied form it is 


LCB =O woes 7. 5C ibis: 
opt Pp 
where we Sais 1s measured in percent L from the forward 
perpendicular. 


Examination of the major series in references (17), (19), 
and (20) which consider LCB as a variable, none adequately 
cover the region in which sailing boats fall. In all cases, 
these series cover block coefficients higher than typical 
yachts. Certain trends may be determined though: 
1. As with the previous factors, there is a peak 
effect at speed length ratio equal to one. 

2. As block coefficient decreases, the penalty for 
non-optimum center of buoyancy position also 
decreases. 


The above series, therefore, place an upper bound on the 


expected penalty. 4A 





Todd (17) indicates for Series 60 that a one percent 
change from the optimum LCB position will increase wave 
making resistance by approximately one percent. Beyond a one 
percent change the penalty increases sharply. 


The assumed factor is 


LCB correction factor = sin (90U,//Z) |62.5-17.5¢,-LcB|"/100 
109 
Frictional Resistance 
The mathematical model for frictional resistance is 
derived from procedures adopted by the Eighth International 
Mowing Tank Conference (21), ITTC in 1957. Traditionally 


mamertyOnal Or Viscous reSistance is defined as: 


a 2 
Re = C,.1/2p Ss, vu, 


S. is the wetted surface of the sailing boat and Ce is the 
coefficient of frictional resistance. 
hiestTTe defined the coefficient of frictional 
resistance as: 
Ce = 0.075/ (log, R - 2) 
n 


where Reynold's Number, Roe is the only parameter affecting 


this coefficient. Again, Reynold's Number is defined as: 
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Ra 7 Up/YH,0 
where Y,0 is the kinematic viscosity of water. 

A suitable definition of the Reynold's Number character- 
@stic length, L, must be established. The -ITTC adopted the 
use of the load waterline length, LWL, as the characteristic 
length. This makes sense for ships with rectangular profiles, 
but sailing boats characteristically have cut away profiles 
and keel and rudder chords which have little dependence on 
the length of the waterline. Davidson (22) realized this 
problem and adopted seventy percent of the waterline length 
as the characteristic length. The adoption of a constant 
fength lacks sensitivity to the effects ciffering profiles 
have on the characteristic length. To make the mathematical 
model sensitive to these variations, the total wetted surface 
of the yacht 1s divided into three areas, each with its own 
characteristic length. 

The three areas are the Sailing boat's canoe body, the 
keel and the rudder. The load waterline length, LWL, is the 
Characteristic length of the canoe body while the mean chord 
of the keel, Che and rudder, Cy are the characteristic 
lengths of the keel and rudder respectively. 

line COetficient Of Eriction for each area is calculated 
using equation 21. The total frictional resistance is the 


Sum of the resistances of each area and is written as: 
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The above coefficients and wetted surfaces are defined as: 


Cech = fairbody coefficient of friction 
Sep = fairbody wetted surface 

Cry = keel coefficient of friction 

Sy = keel wetted surface 

Cee =winigdear coefricient of friction 
S_ = rudder wetted surface 


The ITTC equation for the coefficient of friction is an 
emperical formulation which historically has either slightly 
overpredicted or underpredicted frictional resistance for 
certain forms. From previous experience with a given form, 

a factor may be added to give satisfactory correlation with 

full scale results. When equation 22 is applied to Antiope (3), 
the correlation factor is unnecessary. This would indicate 

that the yacht Antiope has no form resistance or skin roughness. 
This is unfounded and indicates that the ITTC formulation 
Overpredicts the frictional resistance of the yacht Antiope 

by the approximate magnitude of its own form resistance. This 


Meeect 1S discussed further in the next section. 
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Form Resistance 

Form resistance is a combination of effects which show 
little dependence on Froude Number and Reynold's Number in 
the speed range seen by sailing boats. Form resistance, Ree 


is defined as:» 


2 


oi B5)o hs 


fm ~ Cem 17208 


24 


The form resistance coefficient, C is a constant with 


fetyin 
speed and depends only on the boat's geometry, surface 
roughness, and appendages. Form resistance is made up of eddy 
making resistance, roughness effects, and pressure drag due 

to separation. 

The form resistance of a yacht may constitute ten percent 
of the boat's total resistance at low speed and decreases to 
less than two percent as the effects of wave making resistance 
increase. Ea mManieanem Pearl ChmuLs)sai1scuss the difficulty in 
modeling form resistance from analytic and emperical sources. 
The scale effects from small scale model tests have made the 
Measurement of this component difficult. Further, the 
merection Of the ITTC friction coefficient formulation includes 
an arbitrary amount of form resistance which is not quantifi- 
able. The use of this line in predicting Antiope (3) resistance 
did not require the use of a form drag coefficient to provide 


correlation. 
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The present mathematical model does not parameterize 
the form resistance coefficient but allows the user to 
input a value to improve correlation with known results. 
The inclusion of this coefficient allows future parameteriza- 
tion as techniques for its prediction become available. 

Hull Force Formulation 

The total upright resistance, Ru is the superposition of 
wave resistance, frictional resistance, and form resistance 


as defined by equations 12, 23, and 24. 


Ry = RO + Re + Ren 25 
The change in hull resistance with heel and leeway angle must 
be conSidered. For small angles of leeway, yw, the change in 
resistance is considered an effect of lift and is modeled in 
the section on keel forces, 2.1.2. 

In the absence of lift, the change in upright resistance 
due to heel, 9, iS caused by a number of effects. As the 
Sailing boat heels, one would expect not only changes in the 
boat's symetrical lines, but also her form. Depending on the 
yacht form, the effective sailing length may also increase or 
decrease. Further, as the boat heels, the keel moves closer 
to the free surface and in effect decreases the effective 


prismatic coefficient and increases resistance. 
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Examination of a number of towing tank tests conducted 
at M.I.T. show that the increase in resistance at zero side 
force varies depending on the design. This variation ranged 
from an increase in resistance from a minimum of five percent 
to a maximum of twenty-five percent at twenty-five degrees 
heel. The majority of the well-designed models show an 
gaecrease Of approximately ten percent at twenty-five degrees 
heel. Because of the lack of test date available which 
systematically investigates heel effects, the increase in 
resistance with heel is limited to an estimate of what good 
design practice has demonstrated. This effect is emperically 
modeled by dividing the upright resistance by the cosine of 
the heel angle. 

The expression for the hull forces in the Xo Gdirect1on 


™S written as: 


Xy = R/cos > 26 


where Ry 1s the upright resistance defined by equation 25. 


2.1.2 Keel Forces 

The mathematical models for lift forces acting on the 
hull and keel are derived from slender body theory as applied 
to fish-like bodies in references 23, 24, and 25. Expressions 
are defined for the hull and keel side force, Yi induced drag, 


Xi yawing moment, Nii and the heeling moment, Ky 
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Hull and Keel Side Force 

Letcher (3) investigated the possibility of utilizing 
available theories for predicting a sailing boat's side force. 
He compared lifting line theory, slender body theory and an 
adoption of slender-body theory for fish-like bodies to the 
towing tank test results of the yacht Antiope. The theory 
presented by Newman and Wu (24) for fish-like bodies 
demonstrated excellent correlation to Letcher's (3) regression 
analysis of Antiope's test data. An adaptation of this theory 
is utilized for modeling the hydrodynamic lifting forces 
acting on the keel and hull. 

The side force acting on the keel and hull is defined as: 


Dee 


2 
Y, = -C_,cos o(1/20,, ie Uw) 27 


k yw 50 


where the e676 term was found by Letcher (3) to satisfactorily 
gGuantify the effect heel has on the side force coefficient, 
Ey’ for the yacht Antiope. L is the load waterline length. 

A mathematical model for side force coefficient is 
derived from slender-body theory as applied to fish-like forms. 
Milgram (23) gives an expression for the local lift, L(x), on 
a Slender-body as: 


D) 
[b (x) ]V,U XO Coc <x 28 


b 


Syl 





where b(x) is the local span, xX is the horizontal position of 


the maximum span, and x is the horizontal distance from the 


maximum span to the forward perpendicular. Figure 9 shows 
Antiope's profile with the Wea ions ene Xo and x: Vis is 
the component of velocity y direction. 

The total side force predicted by slender-body theory is 


Sb FLO 'O Breb 2 29 


K! 
| 
~ > 


For the numerical computation of the slender-body side force, 
the yacht's profile is tabulated at the five points shown in 
Figure 9 and the profile, b(x), is approximated by straight 
lines between these points. Equation 29 is evaluated 
numerically to give the slender-body lift on the hull and 
keel. 

The velocity he may be approximated for small leeway 


angle by the linear expression: 
ie = UL 20) 


The coefficient of side force, based on slender-body 


“yysb’ 


theory may be written as: 


_ 1 cee 
Cyusb ~ ‘sb/°u.0 & Up¥ ee 


where L is taken as the waterline length, LWL. 
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Newman and Wu (24) developed a procedure for determining 
the effect an axisymetric body of radius ry has on the slender- 
Seretitt. From reference 24, Figure 10 is a plot of 
-L/TpU,V,.b> versus r/o: At rd, equal to zero, the curve 
gives the slender-body lift for a given span. As Ef, 
increases the lift is attenuated by the amount shown in the 
curve. This curve may be defined as a fair body factor, fey: 
The curve is entered with the ratio EO/B. EOVODEaIn Ene stair 
body factor. The slender-body coefficient from equation 31 
is multiplied by the fair body factor to compute the yacht's 


side force coefficient. The expression for the side force 


coefficient is 


Co = Cywsb ‘FED? BZ 
The radius, roe 1s measured as the draft of the canoe body 
at the maximum draft measurement, Do: These two measurements 


are depicted in Figure 9. 

This technique for the solution of side force correlates 
well with the full scale tests of the yacht Antiope (3). 
Antiope has a moderately long keel compared to modern sailboats 
with their fin keels. For a given hydrodynamic span, as the 
chord of a long keel decreases, the lift remains essentially 
unchanged, but at a certain point a further decrease in the 


keel's chord can no longer support the spanwise elliptic 
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loading assumed by slender-body theory. The present model is 
not sensitive to this effect, and though the model 1S adequate 
for modeling changes in draft, it should not be utilized to 
optimize keel area. 

Keel Induced Drag 

Newman in reference 25 develops the equation for induced 
drag on fish-like forms. Figure 11 is from reference 25 and 
gives a curve for (pY,/D;) as a function of the ratio r/D.- 
By generating a spline cubic fit to this curve, induced drag, 


Xp is defined as: 
X, = ¥,w/(pY,/D;) 33 


Keel Yawing Moment 
The yawing moment due to the side force on the keel is 
defined as: 


N 34 


k = Xopk*k 
The center of effort of the hydrodynamic side force, Fook is 
determined from slender-body theory. Milgram (23) gives the 
yawing moment about x in Figure 9 as: 


*h 
Nop = d pol x) ax 35 


O 
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where L(x) is defined by equation 28. Hence, the non-dimen- 


sional coefficient for yawing moment can be written as: 


_ ai Sm 
Causb = Nep/Z°H2° L” ULY 36 


and the distance to the center of hydrodynamic effort from 


the forward perpendicular, x as: 


ela / 


oe = (1.0 -C 


epk )L 27 


ana’ “saints 


Keel He elim Moment 


The heeling moment caused by the side force on the keel 


is 

Ky = ep ew 38 
where Bepk is the distance to the vertical hydrodynamic center 
of effort measured from the waterline. Since the leeway angle, 


Meets smail, cos is taken as unity. 

Slender-body theory tends to overpredict the location of 
the vertical center of effort when compared to towing tank 
data. Letcher (3) reduced Antiope data for the location of 
the center of effort and deduced that the vertical center of 
effort was approximately 1.5 feet below the waterline. From 
this a reasonable estimate for the vertical center of effort, 


Peak’ 1s made. 
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where bo is the maximum draft as defined in Figure 9. 


2.1.3 Rudder Forces 
Mathematical models for the forces acting on the 

rudder are derived from the emperical equation for spade 
rudders presented in reference (26). Expressions are defined 
for the rudder side force, Yi rudder drag, Xi rudder yawing 
moment, ee and rudder heeling moment, K 

Rudder Side Force 

Whicker and Fehlmer (26) conducted a series of wind tunnel 
tests on low and medium aspect ratio spade rudders and pro- 
vided emperical equations which accurately model the forces 
and moments acting on a large family of spade rudders. The 
equations are sensitive to aspect ratio, a, sweep angle, 2%, 
and taper ratio, tL. 

The modeling of rudder forces in relation to a sailing 
boat are complicated by the effect of the downwash from the 
keel. By treating the keel as a low aspect ratio lifting 
Surface, the downwash can be assumed to deflect the flow 
along the centerline of the boat. Although the keel will 
actually deflect the flow slightly less than this, the 
approximation is considered reasonable. The model therefore 


assumes that the angle of attack on the rudder and the rudder's 


deflection from the boat's centerline, a, are identical. The 


effect of heel on rudder performance is considered negligable. 
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Figure 12 is a drawing of a spade rudder including the 
measurements and nomenclature used in determining the rudder's 
performance. The geometrical ratios required by the theory 
are derived from these measurements. The hydrodynamic aspect 
macio, as is defined as: 

ewe 2b. /c, 40 


where b.. is the geometric Span and C. is the mean chord as 


defined in Figure 12. Taper ratio, the is defined as: 


tL. _ Sat rr 41 
where oan 1s the rudder tip chord, and oe Tse tenes rucder, Looe 
chord as defined in Figure 12. 

Whicker (26) defines the coefficient of lift as: 

P| Cc C [Gs 
= L de r 2 
aie Fa) x a ‘Sissy: a 


where Cac is the cross flow drag coefficient, and the rudder 
angle, a 1s measured in degrees. 
An emperical expression from reference (27) for the cross 


flow drag coefficient, C for rudders with faired tips is 


lel 


cdc = 0.08 + O.72t, a3 
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direction to starboard. The z axis desends orthogonally to 
the x-y plane in the direction of the keel. 

The angular differences between the two coordinate 
systems describe the attitude of the boat as it moves through 
the water. The positive direction for these angles are also 
presented in Figure 2. Their direction was selected to pro- 
vide positive values when the boat is in equilibrium under 
normal sailing conditions with a weather helm. The order in 
which these angles are taken 1S important to the final atti- 
tude. The boat is first rotated about the Ze axis through 
the angle ~. The angle wp is identical to the boat's leeway 
angle. The boat is then rotated about the X axis, an angle $¢, 
which is the boat's heel angle. The boat is then rotated 
about the Yo aX1S, an angle 6. 8 is the angle by 
which the boat pitches down at the bow and is also the angle 
of the axis about which the boat heels. Finally the angle Ap 
is the angle the rudder is deflected from the boat's centerline. 

A sailing boat responding to the forces acting upon it 
and conforming hydrostatically to deformations in the free 
surface would heave in the positive Zo direction an amount s. 

The selection of these coordinate systems allows the 
modeling of forces in either coordinate system. When modeling 
velocity related forces, such as wave drag and lift, the 
coordinate system fixed in the water is more attractive 
computationally. On the other hand, hydrostatic forces are 


more easily computed in the coordinate system fixed in the 
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Figure 12- Rudder Geometry 
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aC, /da, is defined (26) as: 


ie, Can9) (27) a 





a eee ee 44 
da. 
16 2g 
(o705)) (COS?) 7 +4 + 1.8] 
COS "9 
where 2 is the sweep angle as defined in Figure 12. 
The rudder's drag coefficient, Car’ 1s defined as: 
C = C + C oy ae+t+ 0.0166 C 2 45 
dr do se laa ; ea 
where e 1s an efficiency factor established at 0.9 (26). 
Ca is the rudder section form drag coefficient and if not 
specified Ca0 is assumed to be 
C30 = 0.0065 46 
mnie 0.0166 sl term 1S an emperical increase in functional 
resistance due to lift recommended by Kerwin (28). This form 


drag coefficient is from the NACA 0015 section (26). The 
interference drag caused by the interaction between the fair- 
body boundary layer and the rudder root chord is not accounted 
for in the present model and should be examined in future 
studies. 

By transforming the rudder's lift and drag coefficients 
from the boat's coordinate system into the system located on 


the horizontal surface of the water, the rudder side force 
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meer ricient is written as: 


Cor = (C, cosy - CarSsinv) cos¢ 47 


Hence, the rudder side force is 


A_U,) 48 
r 


where AL. is the rudder's plan form area. The plan form area 


is defined as: 


Po b C 49 


Rudder Drag 


The rudder drag coefficient is 


Clr = C,,siny 1 Ca pcos 50 


where Cir and Car are defined by equations 42 and 45 


respectively. 
The equation for rudder drag, Xs is 
x =-c_ (5 A_U-) il 


ia xr BNE laa: 
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Rudder Yawing Moment 

To determine the yawing moment caused by rudder forces, 
the location of the center of effort on the rudder, with 
respect to the rudder stock, is calculated. The rudder 
stock 1S assumed to pass through the quarter point on the 
mean chord as shown in Figure 12. 


The moment coefficient about the quarter chord (26) is: 








aC, aC, 
C—,, = [0.25 - (~——), _,] (——._ _ 
mc/4 Te ea ie a 
£ 4B Cac Me: \2 52 
2 a S723 
ie 
where 
3C i I Pe eae 
m meee ‘ ic 
fiGe'c.=0 27 38 
le 4(a +2) 
The normal force coefficient, Caf is 
Ch = C, Cosa, + C5,-COSa, 54 
The location of the center of effort with respect to the 
quarter point on the mean chord is 
C—,,¢ 
> — ee ee 56 
pc/4 C 


nN 
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The distance from the forward perpendicular to the center of 
effort of the rudder, Xopr’ is 


Xcopr a Xoc/4 T Xr 37 
where x,.. is the distance from the forward perpendicular to 
the rudder stock. 


The yawing moment due to rudder forces is 


N = x Y 58 


Rudder Heeling Moment 
From reference 26, the distance to the spanwise center 


me errort from the root chord, Zo! is 





_ (Cy, (473T) cosa, + C, Sina) b, 
z = 59 
ro C 
n 
The distance from the waterline to the center of effort is 
Z, = Z + Z 60 
Gpr pate EEG 
where a is the distance from the waterline to the root 
chord as shown in Figure 12. 
The rudder heeling moment, Ks is 
L i 2 
K- = Zopr“ir ‘2x0 A UL) 61 
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Zoe. 4 Hydrostatic Forces 
A satisfactory model for the hydrostatic righting 
moment 1S presented by Kerwin in reference (5). In modified 


form, Kerwin's (5) model for the hydrostatic righting moment, 


Nee is adopted. 
— 3 4 
N.. = -0.00224[A/(L/100)~] L o (d,+d.o+d,F ) 


+ Zz sing | 62 
2 


where L is the waterline length, LWL, and d d 


17 Go: and d. are 
non-dimensional stability coefficients which depend on hull 
form. The term, Ae Sing accounts for the distance from the 
waterline to the center of gravity. 

Each coefficient has physical significance and may be 
input by the model user or generated from the expressions 
given below. A systematically varied yacht series based on 
the yacht, Standfast 43, is presently being tested jointly 
by M.I.T. and Delft University. The hydrostatic data from 
this series was utilized to generate parameterized expressions 
for each of the coefficients. 

The coefficient, dy, is the height of the metacenter 
above the waterline divided by the load waterline length, LWL. 


An expression for ds which accurately predicts the Standfast 


43 series is: 


4. = 1119 (BWL/L) * 


a 0.335/(L/T.) 63 
[A/ (L/100) ~] 
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The coefficient, d is the non-linear effect on the righting 


oO 
moment. The emperical expression for d., is 
d, = d, [0.0822(BLW/L) - 0.02162] 64 


The coefficient, d., quantifies the effect the wave 
train has on the stability at high Froude Numbers. This 
effect depends on the boat's beam to length ratio, BWL/L, 
and the fineness of her bow and stern. The expression for 
d. is an emperical fit to Standfast 43 data and is only 
dependent on the beam to length ratio. 


5 


[7.8 - 56 ,350(BWL/L-0.28)"' "> (10°) BWL/L > 0.28 


Qu 
tt 
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and 


78% 1077 


Qu 
Il 


SMe = W268 
66 
The location of the center of gravity below the water- 
line, Bas may be specified. If not specified, z_ is assumed 


g 
to be 


Zz. = 0.025L. 67 
| 
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2.2 Aerodynamic Forces and Moments 


The aerodynamic forces are divided into two components, 
the forces acting on the sailplan and its associated rigging, 
and the forces acting on the hull. 

2.2.1 Sail Forces 

The aerodynamic forces are defined in terms of lift 
and drag. Lift acts at right angles to the apparent wind while 
drag acts parallel to the apparent wind. These forces are 
transformed into components fixed in the boat and water 
coordinate system defined in Figures 1 and 2. Expressions 
are derived for the sail driving forces, Xai the sail side 
force, Y oi the yawing moment from the sail driving force, Noy! 
the yawing moment from the sail side force, Nast and the sail 
heeling moment, Ko: 

Wind Triangle 

Computation of the magnitude and direction of the 
apparent wind is necessary for resolution of the aerodynamic 
forces. The convention used in this model was adopted by 


Hazen (6) and defines: 


< 
It 


an Upright apparent wind velocity 


2 2 
¥(V,siny) “- (U,+V, cosy) 68 


and 


68 








(e) 
it 


Upright apparent wind angle 


arctan[(V,siny/(U,+V, cosy) | 69 


where Me is the true wind velocity and y is the true wind 
angle. The trigonometric relationships are graphically depicted 
in Figure 13. This convention differs from the apparent wind 
as measured by instruments on a sailing boat by the magnitude 
of the leeway angle. 

In the case of the sail plan, heel angle changes the 
chord wise direction and magnitude of the apparent wind. 
Since the spanwise component of the wind makes a negligible 
contribution to lift and drag, the apparent wind angle and 


velocity are decreased by the magnitude of the spanwise 


velocity. The apparent wind velocity is 
7 2 2 2 
ey = VV yf (cososing,) + cos ao! 70 
and the apparent wind angle is 


8 = arctan [(cos¢sing,,) /cos8, ] fel 


The resulting lift is normal to the mast. 
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Figure 13 - Apparent Wind Triangle 
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Lift and Drag 


A drawing of a representative sail plan is presented in 
Figure 14 illustrating the measurements utilized by the sail 
force model. 


Lift is defined as: 


eee Ts Fair ve A.) Ue 
where A. is the total geometric area of the sail plan. 
The lift coefficient, Cis! depends on such parameters 
aS sail aspect ratio, boom height, mast diameter, and 
apparent wind angle. The coefficient of lift is also changed 
by crew optimization of trim and camber depending on the wind 
direction and magnitude. The parameterization of the 
coefficient of lift was not done for this model. The model 
user may select the lift coefficient or the model will 
default to 1.25 as a reasonable value. A possible source 
for lift coefficients as a function of sail plan geometries 
is given by Milgram (29). 
The sail drag may be defined as the summation of form 
drag, frictional drag, and induced drag. In coefficient 
form, the expression for sail drag is 
omm 
ve + (5) Le/ 
ls 


i VA O73 


D, = (Cc = Zar a Ss 


1 
Ss fms‘'m * eee ee 





rp 











The friction drag coefficient, C has the assumed (29) value 


ES” 


of 0.04. The form drag coefficient, C based on the cross- 


fms’ 
sectional area of the mast, AW! is a function of the apparent 
wind angle, 8 but to date little data is available on which 
to base a model. For this reason the form drag coefficient 
may be specified by the user or default to a constant equal 


to 0.4. The cross-sectional area of the mast is assumed (29) 


to be 
ae = 0.01(z,, 4. x,,) 74 


where i5 and x, are defined in Figure 14. 
The induced drag of the sail plan is a function of the 
sail plan geometry. Milgram (29) recommends the use of the 


macl1o C oe aS @ measure of sail efficiency. Reference 


dis 
(29) describes a study of systematically varied sail plans 


using lifting surface theory. From this study, an approximate 


Serine cubic fit of C oe was established as a function 


dis 
of the aspect ratio of the sail plan, ass with a correction 


for the difference between the aspect ratio of the main, an 


and the aspect ratio of the fortriangle, a The rae o 


fo 
— é 
Bens’ “1s 1s Gefineda as: 


2. 
@e../"). = f(a) + (a -a,) (0.0625-0.00893 a.) 75 


Or a. 59a 
m fe 


V2 








or 


24. < 
Cais/“is = f(a.) lhe an a. 76 
The aspect ratio of the sail, a., is defined as: 
m2 
= Z5/A, 78 


where Zs is the mast height as shown in Figure 14. The aspect 
ratio of the main is defined as: 

a ao/ Xe 79 
where Z.. and xX, are Shown in Figure 14. The aspect ratio of 


the fortriangle, a is defined as: 


in? 
ap = ne 80 


where Zs and Be are shown in Figure 14. The spline cubic 


Z 


function, f(a.), iomaivonmimeticire 15. The ratio, C /Cy ; 


dis 
is generated for mast head rigs, so the model is capable of 
accepting the ratio designated by the user. The model may be 
used for other geometries by selection of coefficients from 
reference (29). 

The above 1ift and drag coefficients are useful until the 


apparent wind hauls aft and the jib and shrouds interfere with 


the main. At this time the sail plan stalls and a different 
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Figure 14 — Sail Plan Measurements 
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Figure 15 - Sail Efficiency 
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model must be used between the apparent wind angle of 90° and 
180°. Data for systematic variations to sail plans in the 
off-wind condition is not available. The following simple 
form for lift and drag when the apparent wind angle, 8, is 


90° or greater is adopted. 


° a ° 
bog = L.cos (6 - 90°) IO fe! 90 81 


and the off-wind drag is 


D = D. Tae 2 2) 0 


2 ; P 
a A.V. = D.) sin(B-90 ) 


ai 


for g = 90° 22 


If the off-wind lift and drag coefficients are known for a 
given point of sail, they may be specified by the user. When 


data becomes available on reaching the off-wind coefficients 


as a function of the apparent wind, the function may be 


incorporated into the model. 


weewe. wll Forces 
The wind forces on the hull during normal sailing 
conditions may account for only one or two percent in direct 


resistance and ten percent of the side force, but as the wind 


increases and sail is shortened, the hull forces will eventually 


dominate the system. Very few experiments have been conducted 
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investigating the effect of wind on ship forms. Hughs (30) 
conducted tests in water of the topsides of typical tanker, 
liner, and cargo ship profiles. To date no tests have been 
conducted on the sailboat form. Hughs (30) did test the 
above water shape of bare cargo hulls and his emperical 
results are used to approximate these forces on a Sailing 


boat. The magnitude of the aerodynamic hull force, Pur is 


: 2 ee. 2 7 
Fi = 0.001428 We (A, Sin B,t0.3A,cos B,) 7cos (A B) 83 


u 
area of the hull respectively. The angle, A, is the angle 


where A, and A, are the longitudinal and transverse projected 


between the centerline of the boat and direction from which 

the hull force acts. The angle A is taken from the experimental 
results on hulls reported in reference (30). A spline cubic 

was fit through this data and is presented in Figure 16. 

One effect which cannot be ignored when examining aero- 
dynamic forces on the hull is the vertical wind velocity 
profile. The hull of a sailing boat sees considerably less 
wind velocity than the sail plan. Curry (31) measured the 
wind velocity at different heights and plotted a profile. 

By assuming a height of two feet above the water surface as 
representative of an adequate mean for sailing boats, the 
true wind velocity on the hull was estimated from Curry's 


wind profile to be 0.55 times the true wind, V The wind 
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Figure 16 —- Hull Force Angle 
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triangle is calculated with this lower true wind velocity 
and the hull force is calculated with the resulting apparent 
wind velocity and angle. 

Calculation of the yawing moment requires the location 
of the center of pressure of the aerodynamic hull force. The 
distance from the forward perpendicular to the center of 
pressure divided by the load waterline length, was formulated 


from data in reference (30). 
Xepn/ = =e + O20031)] By 84 


The yawing moment, Nue is 


i PLL (x /L)sind 85 


h coh 


Resolution of Aerodynamic Forces and Moments 


The total aerodynamic driving force, Xo is 


Xe = L sing - Dcos8 - Fcosa 86 


and the side force is 


Y. = (L_cosp 1P Dsing)cos¢ Aor, Sank 87 


S h 


We) 





The heeling moment is 


K = Sam) 88 


Z (EEGOSe iebesing) + 2.0 F 
Ss cpys Ss S 


h 
where + ose MomPiometctancceuremchne vertical center of effort 
of the sail side force to the waterline. An expression 
estimating the location of the vertical center of effort, 

, from data in reference (29) is 


LZ 
cpys 


=) .4.5' 2. + Z 89 
aL 


Zcpys £b 


where Zs is the mast height and Z ep is the freeboard measured 
from the waterline to the base of the mast as shown in 
Figure 14. 

The magnitude of the leeway angle is ignored in 
resolving yawing moments. The yawing moment due to aero- 
dynamic side forces, N_, is 


sy 


Noe = Xops (L_cosp + D.sing) cos¢ + N, 90 


where Xops is the distance from the horizontal center of 
effort of the sail plan to the forward perpendicular and is 


defined (29) as: 


Xops dee? Sam + O02 Nx, oa x) 91 
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where Xn is the distance from the forward perpendicular to the 
mast, Xs is the distance from the jib tack to the mast, and 
X is the length of the main boom as shown in Figure 14. 


The yawing moment due to the driving force is 


N = Z X_ sind 92 
Sx cpxs 's 
where 2 eps is the distance from the vertical center of the 
driving force to the center of resistance. Zope is estimated 
by 
Z = 0.48 z. + Zz oo 


cpxs a fb 
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il. GEOMETRY SYNTHESIS 

3.1 Parameter Generation 

The computerized version of the mathematical model has 
the capability to generate dimensions of a satisfactory design 
Beam wiich parametric trade-offs may be made. The user of the 
program may specify as many or as few of the dimensions as he 
desires. This section presents the equations which determine 
the unspecified dimensions. 

The equations do not generate an optimum design. The 
majority of the expressions are estimates by the author of 
the dimensions of a typical ocean racer/cruiser and no other 
justification is made for their formulation. 

If the waterline length is unspecified, the length is 


eet at: 
LWL = 30.0 feet 94 
The default displacement to length ratio, A/(L/100) >, is 
A/(L/100) > = 225 95 
and the prismatic coefficient, oo! is 


C. = 0.53 96 
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The beam to length ratio, BWL/LWL, is defined by a 
linear fit through the mean of a large sample of boats 
presently measured under the International Offshore Rule. 


The beam to length ratio is defined by: 
BWL/LWL = 0.3633 - 0.002018 LWL oy 
By assuming that the base boat has a constant block 
coefficient, Ch! the length to fairbody draft ratio may be 
written as: 


LWL/T_ = 10975 (BWL/LWL) /[A/ (L/100) 7] 98 


As a percentage of the waterline length, the longitudinal 


position of the center of buoyancy, LCB, is 
LCB = 53.0 99 
The fluid properties are established at the following values: 


kinematic viscosity of water: 


a eae ome 100 
2 
density of water: 
= 1.99 Ou 


p 
H,0 
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density of air: 


ian OE UUZ Sic HZ 


The mean chord of the keel, he is estimated by: 


c, = 0.23 LWL 103 


The form resistance and correlation coefficient, Cen! Ss 
assumed to be: 
Cc = 0.0 104 


The estimate of fairbody wetted surface gives a 
satisfactory prediction of the Standfast 43 Series and 
predicts Antiope's fairbody wetted surface to within 1.5%. 
For most modern yacht forms, the expression should give a 
satisfactory estimate for preliminary design purposes. The 


expression for fairbody wetted surface, Sep? is 


= (LWL) *{1.7/ (LWL/T_) +3.5x10~? [A/(L/100) 7] (LWL/T_) } 


{0.0712 [ (BWL/T_)-4.0] + 1.0} ios 


The sail area, Als is established at 


A. = (LWL) * 106 
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The dimensions of the sailplan, as described in Figure 
14, are defined by assuming the aspect ratio of the fortri- 
angle, Aes and the aspect ratio of the main, a, are equal 
to three. The foot of the jib is assumed to be 160% of the 


base of the fortriangle. It follows that the mast height is 
z. = VA /0.4044 107 
at Ss 

and the base of the fortriangle is 


ae = z,/3.9 108 


= 


The main hoist length, eae is assumed to be 


25 = z,/1.1 109 


and the main boom then is 
ILE) 


Zz, = Aa ic 


The horizontal positions of the profile measurements as 


defined in Figure 9 are as follows: 
x = 0.0 1 LAL 


0.1 LWL je2 
85 


m* 
tt 





m* 
It 


5 0.2 LWL i i3 


m* 
il 


3 0.3 LWL 114 


ra 
it 


b 0.55 LWL IES 
The maximum draft, bor is assumed to be the base draft as 
defined by the International Offshore Rule (33) and is 
calculated by: 


bo = 0.146 LWL + 2.0 UALS 


At the other previously defined horizontal positions, the 


draft measurements are 


by = 0.98 bo IIE 
bo = 0.95 LWL/ (LWL/T_) 118 
b. = 0.7 LWL/ (LWL/T.) a, 
bh = 0.0 ZO 


The ratio, ro/bos as defined in Figure 9, is 


ri /B, = Db, LWL/ (LWL/T) lye 
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The sailplan is positioned in the boat so that the 
center of effort of the sailplan, aos! coincides with the 


center of effort of the keel, x as defined in equation 37. 


cpk’ 
The horizontal distance from the forward perpendicular to the 


sail's center of effort, Aoi! is defined as: 





x = x Zo 





| The forward overhang, FOH, is not determined by asthetics, 
_ but coincides with the placement of the jib tack. This point 


is calculated by the expression: 


POH = 0.02 (x. + x_) + X%. - x 125 
J e J cps 


Zep! is derived from the base freeboard 


as defined by the International Offshore Rule (33) with a 
| 
| 


| The freeboard, 


constant added to account for the height of the cabin top 


as shown in Figure 14. The freeboard is defined as: 


Zen = 02057 LWL + 220 | 124 


An estimation of the projected lateral area of the hull, 


Ay: is 


| 

A, = (2.1 LWL + FOH)z,,/2.0 25 
87 

| 
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and the project transverse area of the hull, A is estimated 


t’ 
by the equation: 


A, = 1.1 LWL(BWL/LWL) z 


i ieZ6 


150) 


The wetted surface of the keel, S is estimated by: 


le 
S. = 2.lc [b, - LWL/(LWL/T_)] 127 


k k 


The dimensions of the rudder are determined by the 
following equations. The distance from the forward perpen- 


Gicular to the rudder stock, a is 


x, = 0.95 LWL LZs 
rs 


and the distance from the waterline to the root chord is 

oo OZ LWL/ (LWL/T .) Lins, 
The root chord of the rudder is defined by: 

eee = 0.075 LWL 130 


and by assuming a taper ratio, the tip chord is given by: 


o = 0.667 c 131 
ioe ne 
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The mean chord follows as: 


ee = Cone + chy) /2-0 


An estimate of the span of the rudder is given by the 


expression, 


b,. = 0.85(b, ~ Zz.) 


and the sweep angle, %, is set to zero, 


The rudder area, Ae is 
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IV. PROCEDURE 

4.1 Equilibrium Computation 

To find the equilibrium solution of the forces and 
moments for a given wind direction and speed, equations 1 


through 4 must be solved. The equations are rewritten as: 


R, = X ts xy + X, eS 0 Jig) 7 


Ry = Y. cf vy 4: YY = 0 ii 
R = Ke + Ky 1 Ky + K.. = 0 ILS, 
R, = N aah ees N <0 140 


The solution of these equations require a technique for 
| solving non-linear equations. The technique utilized is 
| Newton-Raphson Iteration described in references (34) and (6). 


R is an "error" vector defined by the components R, 


eechrough Rays 


| R= [Ry R R R 141 


2° Rar Ry] 


where the vector, R, is equal to zero when the sailing boat is 


in force and moment equilibrium. Due to the iterative nature 


of the Newton-Raphson technique, the value of the vector, R, 


is Minimized to an acceptable tolerance. 
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P is defined as the vector of assumed values of the 


four independent variables, 


Pp = [U,, WV, dy ad 


and the first partial derivative 


are defined as: 
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of the component "errors" 
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These partial derivatives are computed by taking first 


meatferences. 


By defining § as a vector of required incremental 


changes to the independent variables, a matrix of linearized 


equations for each "error", R 


mil 12 
po 22 
aon 32 
S41 “42 


13 
Z3 
33 


43 


I 


through R 


is established. 





cm 





Therefore: 


Sin 212 Sag Sn Bias 
| S S S S mie 
(8+ 85s 85, 17 21 822 823 824 2 
Sone 39 oa 234 . <3 
ay sae. Oa -R,| 145 


The new values of the independent variables becomes: 


P, (itl) = P, (i) + 64 146 
P, (itl) = P. (1) 1 65 147 
P,(i+1) = P.(1) 25 53 148 
P, (itl) = P, (2) + Sy 149 


Convergence is tested by examining the magnitude of the 
component "errors", and the vector, P, is updated until the 


desired tolerance is satisfied. 


92 





Vv. RESULTS 

Figures 17 and 18 are performance polar plots of two 
yachts. Figure 17 shows the theoretical boat velocity of the 
yacht Antiope as predicted by the mathematical model. The 
velocity is plotted as a function of the true wind angle, jy, 
and the three wind velocities 5, 7.5, and 10 knots. For this 
Mercdiction thé sail area was assumed to be 400 square feet 
Bnd the center of gravity is 0.61 feet below the waterline. 
The sailplan geometry and placement was generated by the model 
and the keel fixed rudder was modeled by placing the program's 
spade rudder at the trailing edge of the keel and setting the 
wetted surface of the rudder equal to zero. Since the 
correlation of many of the model's forces are derived from 
data of the full scale testing of Antiope (3),-her performance 
polar is included in the results even though the model is not 
designed for accurate prediction of rudder forces on keel 
mounted rudders. 

Figure 18 gives the performance polar of the model 
generated sailing boat for the same three wind velocities. 
The geometry of each boat is fully described in the computer 
printout listed in Section 8.3. The two example output 
listings in Section 8.3 are the source for the plots in 


Figures 17 and 18. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 

Two important criteria should be considered when 
evaluating the capabilities of a system model. The first is 
the verification of the accuracy of the predicted performance, 
and the second iS an assessment of the capabilities and 
limitations a model has for predicting changes in performance 
due to variation in the controlling geometric parameters. 

To determine the accuracy of the model in predicting 
total performance, the performance data of a known yacht is 
required. Unfortunately, the measurement of the performance 
of an actual yacht under sail has never been accomplished. 

To provide a basis for correlation, the four independent 
variables, boat speed, leeway angle, heel angle, and rudder 
angle would have to be meaSured on a full scale yacht along 
with true or apparent wind velocity and direction. This data 
is not available to date. 

Other programs exist for resolving tow tank data and 
assumed sail coefficients into performance polars. These 
programs might be used in the future as a basis for comparison 
if large or full scale towing tank data is used. 

The use of other performance prediction programs for 
correlation is not useful unless another standard exists with 
which to measure their capabilities. If two performance 
programs agree, the agreement may be that both will give the 


same wrong answer. 
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The polar plots in Figures 17 and 18 are considered 
reasonable only because they correlate subjectively with the 
author's own sailing experience. Though each component force 
is correlated against either emperical or theoretical data, 

a means of examining the capability of the marehema tical model 
for the total system should be developed in the future. 

The capability of a mathematical model to respond 
correctly to changes in boat geometry and therefore enable 
parametric tradeoffs requires physical accuracy of the 
individual force models. When the individual force models 
were presented, the limitation and capabilities of each was 
discussed. The limitations and capabilities of each model 


are summerized below. 


Hull Resistance 

Hull resistance is calculated by superposition of 
frictional resistance, wave resistance and form resistance. 

The frictional resistance coefficient is based on the 
ITTC line and is applied to the fairbody, the keel and the 
rudder. This line includes an unspecified amount of form 
resistance. Future model refinement should include a 
frictional coefficient which does not include form resistance 
and a mathematical model that parameterizes form resistance. 
The form resistance coefficient is not parameterized in the 


present model. 
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Wave resistance is estimated by a combination of the 
results of Michell's theory, Taylor Standard Series tests 
and full scale tests of the yacht Antiope. The resulting 
semi-emperical mathematical model for wave resistance needs 
to be compared to the results of a series of full scale yacht 
wave resistance tests. 

The mathematical model for ‘hull resistance is considered 
adequate for anadliotanang calm water resistance. To predict 
a boat's performance in rough water, a mathematical model 
predicting the added resistance due to waves should be 
developed. 

The change in resistance due to heel angle is also not 
parameterized and was modeled as an average of what good 
designs can expect. The investigation and incorporation of 
heel effects would be worthwhile. 

The present resistance model is considered dependable 
for parametric variation in displacement to length ratio 


and wetted surface. 


Keel Forces 

The hydrodynamic side force on the hull and keel are 
derived from slender-body theory modified for fish-like forms 
(24). This model provides accurate prediction of forms fixed 
with moderate keels such as found on the yacht Antiope. A 
method for determining decrease in lift as the keel chord 


decreases at a given span should be developed. The present 
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method adequately models changes in performance due to draft 
but is insensitive to reduction in keel area due to chord 
length changes. 


The interference effects between the hull boundary 


player and the root chord of the keel should be the subject 


of further investigation. 


Rudder Forces 
The rudder forces are modeled from the free stream tests 


of low and medium aspect ratio lifting surfaces (26). The 


| model is adequate for yachts with spade rudders, but many 


yachts are being designed with skeg rudders. A model should 


be developed for rudder forces as a function of the percentage 


flap on the skeg rudder. The fairbody boundary layer inter- 
ference effects on the rudder have been ignored in the 
present model. 

The mathematical model is expected to satisfactorily 


model changes in the spade rudder geometry. 


pydrostatic Forces 


Hydrostatic forces are modeled uSing the concept of the 
metacentric height and are adjusted for Froude Number effects. 
The model is considered sensitive to geometric changes though 
further verification of the Froude Number effects is desirable. 

To accurately model the righting moment on smaller yachts, 
a function representing the effect of crew weight on righting 


moment should be added. 
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Ball Forces 

The sail force mathematical model is based on lifting 
surface theory presented in reference (29). The model is 
adequate for close hauled prediction but is not parameterized 
for reaching and running due to lack of data. Work is needed 
to determine force coefficients as a function of the apparent 
wind angle and sail plan geometry for off wind conditions. 
Further work needs to be accomplished to parameterize the 
form drag coefficient as a function of rig geometry and 
apparent wind angle. 

The utilization of the model should be limited to light 
and moderate wind conditions until the optimization variables 
for reefing and flattening are included. 

The effect of the vertical wind gradient on heel has 
not been included and should be added in the future. When 
added resistance in waves is investigated, it will also be 
important to study the effect large pitching motions have on 
sail forces and their relationship to the added resistance 
in waves. | 

The limitations of the sail force model in its present 
form preclude the accurate prediction in off wind and heavy 


air conditions. 


Beredynamic Hull Forces 


The wind forces acting on the hull are modeled from 
data on cargo ship hulls (30). Though the model is derived 
from a different hull form, the prediction of these forces 
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is considered reasonable. Tests on yacht hulls would help 
improve the correlation for this mathematical model, but 


the effect on the final prediction should be minor. 


In Conclusion 

The present form of the mathematical models is useful 
for preliminary design in determining the quantitative effect 
changes in the major geometrical parameters have on a boat's 
performance. These parameters include displacement to 
length ratio, wetted surface, beam to length ratio, total 
draft, rudder geometry, position of the center of gravity, 
sail plan geometry for on wind conditions, and sail area. 

The development of the model has been useful in cata- 
loguing available theoretical and emperical sources for 
modeling forces acting on a sailing boat and has made clear 
the necessary improvements required in the science of sailing 
boats to improve prediction capability. For the naval 
architect, this first generation system model for sailing 
boat performance may be used as a file to which improvements 
may be made as a greater understanding of the science 


develops. 
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merit. APPENDICES 
@reecomputer Input 
The required data cards are divided into two data 


sets described below: 


Card Set Number One 

Card set number one consists of the input points 
necessary to generate the spine cubic curves required by the 
force models. There are twenty-four cards punched in FORMAT 
(8F10.0). This data set is an integral part of the program 
and does not change from run to run. Table 1 is a listing 


of the required data cards. 


Card Set Number Two 

Card set number two 1S repeated for each boat 
evaluated by the mathematical model. The cards control 
the program options and describe the geometry of a given 
sailing boat. A listing of the card set used to describe 
the yacht Antiope is given in Table 2. Table 3 is an 
example of the required data set where all parameters 
are defaulted. 

Card Set #2 

card #1 IFAM, NGAM, NVT; FORMAT (8110); no default. 
sol 10 If IFAM = 1 forces and moments are printed. 


If IFAM = O no forces and moments are 


printed. 
eel li - 20 NGAM = number of true wind angles desired. 
col 21 - 30 NVT = number of true wind velocities desired. 


OS 











Sard #2 


mol lil - 13 
Bol 31 

eol 53 

card #3 

eo. lil - 17 
Sol 3l - 37 
card #4 

col 21 - 30 
col 31 - 40 
col 41 - 50 
eol 51 —- 60 
card #5 

or more 

col 1 - 80 
card #6 

or more 

col 1 - 80 


Pie Selenite Ns PORMAT (1OX, 13, 17X, I1, 
Uses) > no default. 


ITMAX = maximum number of iterations 
allowed per solution. 


If IPRINT = 1, the independent variables, 
apparent wind data, and forces and moments 
will be printed out each iteration. 


If IPRINT = 0, only the convergent or 
final iteration will be printed. 


N = 4, the number of degrees of freedom 


Epo, EPSZebORMAT (10X, E7.1, 13x, 
E/.1); no default. 


EPS1 = required tolerance in force 
eyo bial Urtlepentibigi< 


EPS2 = step size taken for computing 
first differences. 


CXC CL) ata eer ORMAT (20X, 5SF10.3); 
no default. 


XOLD(1) = initial assumed value of boat 
velocity, (feet per second) 


XOLD(2) = initial assumed heel angle, 
(degrees) 


XOLD(3) = initial assumed leeway angle, 
(degrees) 


XOLD(4) = initial assumed rudder angle, 
(degrees) 


Cowl) te—weniel) > FORMAT (8F10.5); 
no default. 


KVT(L) = desired true wind velocities 
(knots) 


OGAMMA(I), I = 1,NGMA); FORMAT (8F10.5), 
no default 


DGAMMA (I) = desired true wind angles 
(degrees) 
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mara #7 


oo 1 — 10 
eo 11 - 20 
eel 21 - 30 
mol 31 - 40 
Sot 41 - 50 
card #8 

eo. 1 - 14 
card #9 

esol 1 = 10 
Bel 1! - 20 
esol 21 - 30 
col 31 - 40 
pot 41) - 50 
gor 51 - 60 
col 61 - 70 
card #10 
col 1 — 10 
Sor ii - 20 
Bot 21 — 30 
sol 31 - 40 
sol 41 ~ 50 


C(2), 
(8F10. 


G2) 
C(7) 
G8) 
Ci) 
€ (20) 
Grell), 
exert} 


(C(I), 


Ciwomeeus), €C(9)- C{10)> FORMAT 
5); default = -1.0 


waterline length (feet 
displacement to length ratio 
prismatic coefficient 
lemgeh to draft ratio 
= distance from FP to LCB (%LWL) 
FORMAT (E14.7), default = -1.0 
= kinematic viscosity of water 


iv—eelee, Oo) > FORMAT (8F10.5); 


default = -1.0 


CxX12) 
GES) 
C(14) 
GES) 
GEG) 
OIE), 
eC) 


GO 
(sila 8) 6 


density of water 

= mean keel chord (feet) 

= mean rudder chord (feet) 

= form and correlation coefficient 
= fairbody wetted surface 


= keel wetted surface 


rudder wetted surface 


eleven s) ec (5) ,-C(6)> FORMAT 
5); default = -1.0 


C(1) = beam to length ratio 


C(3) = sail area 
C(4) = heeling moment arm 
C(5) = pitching moment arm 


C(6) = distance from FP to rudder stock 


HAG) 7 





card #11 


Zot | - 10 
mee ll - 20 
oe 21 - 30 
mort 31 - 40 
mol 41 - 50 
mee 51 —- 60 
men 61 - 70 
fol Jl - 80 
card #12 
oor 1 - 10 
mee ii - 20 
Bor 21 - 30 
fol 31 - 40 
card #13 
fol 1 - 10 
pov it — 20 
Sol 21 - 30 
eel 31 - 40 
fot 41 - 50 
mea 51 - 60 
pol 61 - 70 
cot 7l — 80 


(GL). Tr 
default 


C(19) = 
20) = 
€(21)" = 
G22) a= 
C( ZR) = 
C(24) = 
es). — 
C(26) = 


(GAGE)... 5 
default 


CZ 7) 


EZ 3) 


eZ 9) 


Ce 0) 


(ey LS 


default 
Csi) = 
Gie2) = 
Ges) = 
C(34) = 
E35) = 
C(36) = 
C(37) = 


C30) 5 


iheaoue FORMAT (SF10.5) > 
Se) 


close hauled lift coefficient 
rig form drag coefficient 
cross sectional area of rig 
Sail frictional coefficient 
induced drag coef/lift see 
lateral area of the hull 
transverse area of the hull 
auseaneces=ron Ch sall to FP 


Zoe eer ORMAT (SF 10.5); 
alee 8, 


off wind drag coefficient 
off wind lift coefficient 
Gistanccmmen Chahuli to FP 
density of air 


Soom EORMAT (8F10.5)> 
= -1.0 
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card #14 


Zot Lt - LO 
mot ll - 20 
card #15 
fol Lt - 10 
fot 11 - 20 
Zon 21 - 30 
card #16 
fot 1 - 10 
fot 11 —- 20 
mot 21 - 30 
Zol 31 - 40 
sol 41 - 50 
Bol >i ~ 60 
Bol 61 - 70 
col 71 - 80 
eard #17 
sol 1 - 10 
mor ll —- 20 
mot 21 — 30 
eol 31 - 40 
fol 41 - 50 


Cee e (40) ePORMAT (S8F10.5); default=-1.0 


39) b 


s) 
Dy 


C (40) 


(C(I), I = 41,43); FORMAT (8F10.0) ; 
default = -1.0 


C(41) = draft of the canoe body/total draft 


C (42) 


hydrodynamic lift coefficient 


G43) 


WX lift/drag 


(enh) pets) FORMAT (8F10.5) > 
default = -1.0 


C(45) = distance from LWL to rudder 
root chord 


C(46) = rudder root chord 
G47) = rudder tip chord 
C(48) = rudder span 

C(49) = rudder sweep angle 
C(50) = rudder aspect ratio 
C(51) = rudder mean chord 


G52) 5 — ruader 9C/ 90. 


eae te —e obo) )-C(61)> FORMAT (8F10.5) ; 
default = -1.0 


Gi5s))—= ruader cross flow drag coefficient 
Gie4— rudder aC _/9C, 

C(55) = rudder form drag coefficient 

C(56) = rudder area 

C(61) = distance from CE keel to LWL 
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card #18 (G(iies ft = 57,60); FORMAT (8F10.5) > 
default = -1.0 


col 1 - 10 C(57) = metacentric height/LWL 

Bot li - 20 C(58) = Froude Number righting moment 
coefficient 

mol 21 - 30 C(59) = non-linear hydrostatic coefficient 

Bol 31 - 40 C(60) = distance from LWL to CG 

card #19 (Geb i —s62,69)3 FORMAT (8F10.5); 
default = -1.0 

mol L - 10 C(62) = mast height 

col il - 20 C(63) = base of the fortriangle 


mole2zl — 30 C(64) = main hoist 
col 31 - 40 C(65) = main boom length 
col 41 - 50 C(66) = aspect ratio of the fortriangle 
col 51 - 60 C(67) = aspect ratio of the main 
Sol 61 - 71 C(68) = aspect ratio of the sail area 
sot 71 — 80 C(69) = distance from mast aft to sail CE 
card #20 Ce OvemCivi) merORMAT (8F10.5); default=-1.0 
eo! 1 - 10 C(70) = forward overhang 
wei il — 20 C(71) = freeboard 
card #21 KEND, FORMAT (I1), no default 
sol lL If KEND = 1, indicates this is final 
data set 


If KEND = 0, indicates other data sets 
GOmrEO. LOW 
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8.2 Computer Listing 


The main program and all necessary subroutines 
utilized to compute the mathematical model are listed 


in this section. 
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((VLLG9G) SODxOVAT- (WLEAG) NISeLIIT) = (Lt) A 


SPOWIOd IONTATYG OAV TLINdWHOD 


ANNILNOD 
(€) DxZ0% (27) D=ZOWUG 

(€) DeOWe (22) D=9VECY 

(€) D*O% (LZ) D=9DVUG 

0Z OL O9 
(LG°L-ZVELIG) NISx (ZOVHa— (EF) DeZOu7Z OL) +ZOVUC=ZOVIG 


(LG°L-VLLIGY) NIS« (OVNaY- (€) Oe OVeZ°L) +OVEGV=9VHAaY 


—-- 


0¢ 


Lt 


UO O OoN YU 
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myo <) 





ang 


NUNLaY 
(te) A=(h’ Lb) AW 

(LL) A=(E“LLD AV 

((Z) GTOXWINIS« (ZL) AV& (GS) 9=(Z LL) AY 
((Z)GTOX)NISs (L4L) AW (IG) D= CL ‘ELD AY 
((Z)GTOX) NISx (QL) Ax (GJO=CLL)A 


ADNYLSISINY dO USLNAD OL IIVS AZOWONA ADNVISIAG =(SG)9 
(tL) % “INGWOW ONIMVA TIVS ALNdAWOD 


(OL) A=(h’OL) AV 

(OL) A= (EOL) AY 

HWA+ (9Z) De (( (L) VHdTV) NIS#Hd-(7°S) AW) =(Z’ OL) AY 
HWAV4(9Z) Oe C((Z) VHa TW) NISenaAV- (LS) AW) = (LOL) AV 
HWA+ (97) De (((L) VHATY) NTSxHI-(S) A) = (OL) A 


Sd940d ACIS OWVA OL ANG ILNAWOW ONIMVA ALNdHOD 


(g)A=(h'9) AW 
(9) A=(€'8) AV 


(tr) Del 


((ZVLELAG) NTIS «Z9OVNaG+s (ZVLLIG) SODeZLAI1) = (7° 8) AV 


(fh) Del 


OU2UDY 


OO YO 
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Zak (0°7-( (1) d99H) OLDOTY)/SL0 °O= (I) aua2 

Cx (0 °2-( (1) TY) OL DOTY) /S20 °O= (1) TaD 

(GL) 54Ze% (O°¢- ( (1) da34) OLDOTY) JGL0 °O= (1) daad 
O°L=(T)auga (o°0°sST° (1) Guay) aT 

OT b= (2) tae (0°07 8a) oa ea 

O°t=(I)aaga (0°O°S1° (1) aGagu) al 

(LL) O/4 (1) Sax (ht) o=(1) auay 

(LL) O/ (1) SA« (€L) 9=(1) THaY 

(LL) D/ (I) SAx (7) D=(1) Gagy 


SONVESTISANY TVNOILOIYd FYLNdWOD 


((1) dod O41 

(I) H1d94 (1) d9dd940°L) *(L) O% (4) Daf axe (0° O00L/(Z) 3) « (1) La= (4) dAVG 
O°OOL ZExe (( (OL) D- (8) D%S °LL-S°729) SAV) & ( (I) SIS*LG6°L) NIS=(I) dd1dd 
((6)9-0°9L) «(r) €4a-= (1) HTa9 

L°txx ((O°OOL* ((1) dOdd- (g)5)) sav) * (1) Za=(1) ddad 

Zt eeoG 


JINVLSISHY JAWM ALNdIWOD 


(L)aTOXV=(Z2) SA 

(L}aTOX=(L) SA 

(dav ’dOdD°HIS4NIX‘4Z4LL) NATWAR TIVO 
(CAVEA*MIS“NIX“Z“LL) WATWAZT TIVO 

(Tau Ca‘UIS‘NIX‘’Z‘LL) NYATYAR TIVO 
(LaLa HIS“ NIX’ CLE) SQUAD ao 
G°Oax(Z) 9/689 °LZ (L) ATOXV¥=(2) HIS 

G°O«« (7) 9/7689°L/(L) aTOX=(L) aTS 

(7) aTHYaa’ (2) auad’ (Z) THAD’ (7) GaqD NOISNAWIA 

(7) SA’ (Z) 08TN NOISNANIGC 

(Z)tngu ’(Z) aagyu’ (7) aavaua’ (Z) ga74d904 (Z) Hdd’ (Z)dDID NOISNAWIA 
(7) d0d9’% (Z) €47% (2) Za’ (Z) La’ (Z) UTS NOISNAWIG 

(LL) NIX’ (hh) adv’ (hh) ea’ (hh) 7904 (ho) LAY NOISNAWIG 
(h“GL) AW’ (GL) 2% (O0L) 94 (4h) ATOXV’ (h) aGTOX NOISNAWIAG 


(NT X’dHOV’EAV ‘TAN LAY AW‘ A4OD4aTOXY‘“GTOX) TINH ANILAONGNAS 


OU YO 
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ang 


NYUNLAY 
(2) R= (Gee 
(2) h=(6 -o) A 
(((z)aToxv)Ssoo/ ((L) oTuaas (L) tAVMG)) Sav-=(7'Z) AV 
(((Z) d10x) Soa/((z) DTH¥aas (7) TAWNa) )SAV-=(L‘“Z) AY 


SAONOT AGZLINAWONY ALNdkOD 


(((z) adTOX) SOD/ ((L) DTHaas (L) gawMa) )saqv-=(Z) a 
AINNILNOD Z 
((8L) Det 


(I) q¥do+ (LL) Dx (1) THAD4 (OL) Dw (I) PAID) Zee (IT) SAKO°C/ (ZL) D=(T) OTN 


UOUUY 
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Zxa¥ x (0S) D4 (ES) 9«S °O~ (ZS) ae ( (nS) D-Sz°0) =HD 
(VVINITS*hdo+ (¥¥Y)SOD*xtfTID=hNo 
(¥)NIS%d9+(¥)SOD«TID=N] 


(CL) A “USdadNY WONT LNAWOW ONIMWA ALNAWNOD 


Ox ((€X) SODehaa+ (EX) NISxhT9)-=(h’ th) AY 
Ox ((EXX) SOD*dD+4+ (E XX) NIS«T9) -=(e% th) AY 
(nh) A=(7'% tH) AV 

OVeOVG-=(b’ h) AV 

OxIVI-=(h) Z 

(€X) SODedD4+ (EX) NIS&xTDO=OVd 


(mh) A “OVE UITANHY ALNanod 


((Z7) dTOX) SOD«O x ( (EX) NISehdd- (€X)SOD*nITD) -=(h’Z) AY 
((Z) dTOX) SODeOe ( (EXX) NISad90- (€XX)SOD*%T3) -=(€°L) AV 
((Z) GTOXY) SODKO«OVI-=(Z°L) AV 

((Z) GITOX) SOD*OV*OVd-=(L OL) AY 

((Z) @TOX) SODaO «DVI -= (2) K 

(€X) NISexd9-(€X) SODxTD=OVd 


(L)A*TDNOd AGIS WIACANH ALNdwod 


(9G) DaZax (L) GIOXVx (ZL) De G* 0=0¥ 

(9G) De Zxe (L) GTOX x (ZL) DS °0=0 

CeaTTIRIOLO °04 (105) 9%L78°7)/7x«NTD+ (GS) D=had 
Ce*1T9%991L0°04 ((0S) DeL78 °Z) /Zx%xT94 (GS) D=d9 
Ze%V x (0S) 9/7 (ES) 24 (h) aTOXVx (ZS) D=hT9 

Tx*WVx (0S) D7 (ES) D4 (h) ATOX« (ZS) D=19 

CC 7 Ae) dior oy 

£€°2S/ (€) dIOX=EX 

E°LG/ (hh) GTOXV=V¥ 

€°£S/ (hn) GTOX=V 

(h°St) AWS (GL) A*% (COOL) D4 (h) GTOXY’ (fh) ATOX NOISNANIG 


(AV ’°A*D ‘aATOXY’ATOX) HaaqnY ANILNOUGNS 


UU DO 


OO Nv 
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GNG 


NUNLaY 
(Sdd4+ (St) 9) sOxthTo= (h“SGL)AY 

(Se) A= (e* SpA 

| (SL) A=(Z°SGL) AW 

(SdD+ (Gt) D) xOVxTO= (L“SL) AY 

(SdD4 (Gh) 9D) eOxTO=(GL) A 

(8th) DethNO/ ((WY) NTS %hGd+ (YY) SODxth1D«hh7H°O) =hSdd 
(8t) De NO/ ((Y) NISeI94 (¥) SOD*%TO«thnZh°O) =Sdd 


YAdGNY WON TNEWOW ONTITAIGH ALNdWOD 


(hdd+ (9) 9) x (nL) AV= (HEL) AV 
(79+ (9) O) x (ELI AW=(E°EL) AV 
(dd+ (9) 9) x (Z4L) AV=(Z°EL) AV 
(do+ (9) OD) we (LL) AW=(CL ELD AY 

(49+ (9) 9) x (LJ AH(EL)A 

(LS) Ox NO/thWO-=hdd 

LOOO°O=HND (LOOO°O"°ELI*HND) aT 
(LS) D&k&NO/WO-=dd 

LOOO°O=ND (LOOO°O °LI‘ND) AI 


J wxVV x (0G) 9/ (ES) 9%G°O- (7S) De ( (HG) D-G7°0) =hWa 
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LNIWOW ONITIAH ALAdWOD 


(TL) A= CH°TLI AY 

(€°9) AV%(62Z) D= (€ °ZL) AY 
(779) AV % (67) D=(Z°ZL) AV 
(1°99) AWx (67) D= (LTE) AY 
(9) X*« (62) D= (74) A 


(CL)A “ADHOd AGIS THAW OL ANG INIWOW ALNaAWOD 


(€)k=(h‘€) AV 

Tux ((Z) GTOX) SOD/ (En) D/E°LG/ (E) CTIOXWeE TX-=(€°F) AV 
Tx ((Z) TTOXY) SODS (EH) OZE °LG/(E) ATOX*ZTX-=(Z°E) AV 
- xx ((Z) GTOX) SOD/ (Eh) OVE LSS (E) GTOXxLTX-= (LE) AV 
7% ( (7) GT0X)SOD/S (Eth) OZE °LGZ (€) GTOXx1X-=(E) XK 


(c) A *° OVA IFIN GIDNINI ALNAWOD 


(9) X=(n‘9) AV 
SOO CaS ey 
CTX-=(7°9) AV 
LTX-=(L’°9) AV 

TX-=(9) i 


(o)k “TONoOa ACIS TARY BFLNAwWOoOd 


Cae (Z)OKeE LG (CE) ATOXV xO (Zh) D=E 1X 

Cx% (Z)DeE °LG/(E) ATOX*ZO% (7H) D=7TX 

Cx (Z) Oe €°LG/(E) GTOX xt Ox (Zh) D=L7X 

Txx (7) De€ °LG/ (€) GTOX*O« (Zt) D=1X 

0° 7/Zx x (((Z) dIOX¥) SOD» (L) @TOX) x (ZL) D=ZO 

0°7/Ze« (((Z)0T0X) SOD« (L) dTIOXW) « (ZL) D=L 0 

0° C4 Zax (((Z) 0T0X) Sod* (L) GTOX) x (Zt) D=0 

(nh?°SL) AV’ (SL) A% (001) 9% (fh) aTOXY ’ (hn) aTOX NOISNAWIG 


(AV°AR’O“GIOXW“QTOX) THAW ANILAONNS 


ees) CESS. SBOE: 
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Gna 


NYOLaY 

(hl) A= (h’ HL) AY 
(o°L—-) « (07) adTOX) Soa/ (649) AV (19) O= (E “hL) AY 
(oO°L—) x ((Z)aTOXW)SODZ (7 °9) AV (L9) D= (27 HL) AY 
(O°t—-) x ((Z)GTOX)SOOS (L°9) AV (L9) D= CL “LD AV 


(ocLr—-) « €(Z) a10X) SOD/ (9) Ax (19) D= (th) 2 
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ana 


NUNLAY 
(6) K=(’6) AV 

(6) k=(€°6) AW 

((Z) GTOXW) NISx (09) D414 

((Z) CTOXV & (6S) D4Nde (8S) 04 (L5)9) x (7) CTOXW* OVI=(7'6) AV 
((Z) G1OX) NISx (09) D¢L 

((Z) GTOX x (65) D4 Nd x (86) 94 (25) 5) « (7) GTIOXxDVa= (E46) AY 


SHOYOd AFTLNAWONW ALANdWOO 


((Z) ATOX) NISx (09) D4L 

((Z) GTOX « (696) D4+Nd% (85) 94 (25) 9) « (Z) GTOXx DVa= (6) K 
((Z) De 72°ZE) LHYOS/ (L) ATOXV=LNGA 
((Z)D«Z2°ZE) DIOS (L) GTIOX=Nd 

O°tex (Z) Ox (L) Det2 700 °0-=9Vd 

FAT = (Z)d 

TMT / Tmad =(1)9 


LNGWOW ONILHSINY OLLVLSOUGAH RLINdWOD 


(n’St) AW’ (GL) A’ (COOL) OD’ (fh) aTOXW’ (hH) aTOX NOISNAWIAG 
(AV ‘(RO “ATOXY’ATOX) NONLY SNILQONGNS 


CUS EN 


OOUU0U0 0 
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8.3 Example Output 


Two example computer outputs are listed in this section. 
The first is the yacht Antiope and the second is the model 
generated sailing boat which is produced if all geometric 


parameters are defaulted. 


lay 





ITsaxr = 16 
IPRINT = 6 
u = a 
BPS1 = 1. OD- 10 
EPS2 =, 1. OZ-02 
XOLD(1)..-KXOLD( 4) 
1.000000E 01 0.000000E-01 0.000000E-01 0. OOO000E-01 


SOLUTIONS ARE COMPUTED FOR THE FOLLOWING TRUE WIMD SPEEDS (KNOTS) 


5.00000 7.50000 10.00000 ; 
SOLUTIONS ARE COMPUTED FOR THE FOLLOWING TRUE WIND ANGLES FOR EACH WIHD SPEED 


180.00000 160.00000 140.00000 120.00000 100 .00000 80.00000 60.00000 45.00000 35.000006 
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LWL= 24. 30000 

DISPLACESENT TO LENGTH= 180.39990 
PRISMATIC= 0.53000 

L /a= 12. 15000 


LCB= 53.50000 
KINEMATIC VISCOSITY= 0.1279100E-08 


DENSITY OF WATER= 1.99000 

MEAW KEEL CHORD= 7.12500 

BEAN RUDDER CHORD= 0.00000 

ROUGHNESS AND CORRELATION FACTOR= 0.00000 


PAIRBODY WETTED SURPAC E= 118.61000 


KEPL WETTED SOURPACE= 41, 38200 
RODDER WETTED SURPACE= 0.00000 
TOTAL WETTED SORFACE= 159.99200 
BYL /LWL= 0.23900 


SAIL AREA=  400.00000 

HEELING MOMENT ARM= 16.27361 

PITCHING MOMENT ABRA= 17. 84613 

DISTANCE PRON PP TO RUDDER STOCK= 16.00000 
CLOS2 HAULED LIFT COEF= 1.25000 

RIG PO&8S DRAG COEFP= 0.40000 

CROSS SECTIONAL AREA OF RBIG= 19.60001 

SAIL FRICTIONAL COEF= 0.04000 

INDUCED DRAG COEF / LIPT COEP**2= 0. 13030 
LATERAL AREA OF HOLL= 100.00000 

TRANSVERSE AREA OF AULL= 27.00000 

DISTANCE BETWEEN CE SAIL AND PP= 11.50000 


OPP WIND DRAG COEFP= 0.00000 


Jasye: 





OPP WIND LIFT COEP= 0.00000 

DISTANCE BETWEEN CE HULL AND PP= 11.3576 

DENSITY OF AIR= 0.00238 

HORIZONTAL DISTANCE PROM MAXISUS DRAFT OF PROFILE MEASUREMENTS 
0.00000 5.75000 8.51000 14.67000 16.91000 

PROPILE AEASUREMESTS 


4.65000 “4.62000 2.12000 0.78000 0.00000 


DRAPT OF CANOE BODY / TOTAL DRAFPT= 0. 35000 
S 
HYDRODYNAMIC LIFT COEP=. 0.09805 
LAMDA X LIPT / DBAG= 2.06498 
DISTANCE FROM LWL TO VERTICAL HYDRO CE= 1.598495 


*** RUDDER SECTION *## 


DISTANCE PROM L¥L TO ROOT CHORD= 0.80000 
ROOT CHORD= 1.82250 

TIP CHORD= 1.21561 

RUDDER SPAN= 3.61250 

SWEEP ANGLE= 0.00000 

ASPECT RATIO= 8.75625 

MEAN CHORD= 1.51905 

DCL / DALPHA= 0.06795 

CROSS PLOW DRAG COEPF= 0.56024 
DCH / DCL= 0.21407 

PORN DRAG COEF= 0. 00650 
RUDDZR AREA= 5.48758 


=e RIGHTING MOMENT SECTION #22 


METACENTRIC HEIGHT / LYL= 0.05711 


PROUDE NUSBER RIGHTING MOMENT COEFP= 0. 78000 00E-08 
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NOM LINEAR HYDROSTATIC COEP= 


0. 112748S1E-03 


DISTANCE BELOW LWL TO cG= 0.60750 

** SAIL SECTION *## 

I= =. 31, 45027 Jz 10. 48342 

P= =. 28. 59118 Re 9.53038 

AR POR TaI= 3.00000 ° AR MAIN= 3.00000 
AR SAIL PLAN= 2247280 

CE SAIL PLAN APT BAST= 0.90028 


PORWARD OVERHANG= -~0.61630 


FREE BOARD= 2.75000 ° 
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KYT= 5.00000 KB8OTS DGABSA= 180.00000 


KYA= 2.48617 KNOTS DBETTA= 179. 99980 


SUCCESSFUL CONVERGENCE ° 


ITER = 5 


BOAT SPEED= 2.51383 KNOTS 


LEEWAY ANGLE= -0.00390 DEGREES 


a 


x 
eg ebhLeouee ce @& Se seen Feee2 ee Feeweeaeneaae eS Gee eet aneroe ee ereo om 


HEEL ARGLE= 


RUDDER ANGLE= 


DEGREES 


DEGREES 


0.00128 DEGREES 


0.01191DEGREES 


KYT= 5.00000 KNOTS DGABSA= 160.00000 DEGREES 


KYA= 2.60820 KHOTS DBETTA= 139.03030 DEGREES 


SUCCESSPUL CONVERGENCE 
ITER = 2 
BOAT SPEED= 


LEEWAY ANGLE= 0.02088 DEGREES 


NO CONVERGENCE 


KYT= 5.00000 K8OTS DGASHAA= 140.00000 


KVA= 3.28689 KSOTS DBETTA= 
SUCCESSPUL CONVERGENCE 


ITER = 11 


BOAT SPEED= 


LEEWAY ANGLE= 0.148458 DEGREES 


2.72703 KWOTS BEEL ANGLE= 


RUDDER ASGLE= 


98.17028 


3.36880 K¥oOTS HEEL AYGLE= 


RUDDER ANGLE= 


0.02578 DEGREES 


0.01149 DEGREES 


DEGREES 


DEGREES 


0.40965 DEGREES 


0.00311DEGREES 


KYT= 5.00000 KNOTS DGASMA= 120.00000 DEGREES 


KVA= 4.67690 KSOTS DBETTA= 


SUCCESSFUL CONVERGENCE 


56 


67.79747 DEGREES 





ITER = q 


BOAT SPEED= 4.26728 KNOTS HEEL ANGLE= 2.74890 DEGREES 
LEEWAY A¥GLE= 0.41392 DEGREES : RUDDER ANGLE= 0. 19096DEGREFS 
KVT= 5.00000 KNOTS DGAASA= 100.00000 DEGREES 

KVA= 6.28702 KNOTS DBETTA= 51.55513 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER =z § . 
BOAT SPEED= 4.77718 KNOTS HEEL ANGLE= 6.77472 DEGREES 
LEEWAY ANGLE= 0.66322 DEGREES RUDDER ANGLE= 0. 66225DEGEEPS 
KVT= §.00000 KNOTS DGASMA= 80.00000 DEGREES 
KVA= 7.53630 KNOTS DBETTA= 40.79662 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = 3 


BOAT SPEED= 8.83670 KNOTS HEEL ABGLE= 10.99982 DEGREES 
LEE@AY ASGLE= 0.97314 DEGREES RUDDER ANGLE= 1. 19417 DEGREES 
KVT= 5.00000 K¥OTS DGAHMA= 60.00000 DEGREES 

KVA= 8.15836 KNOTS DBETTA> 32.05679 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = 3 


BOAT SPEED= 4.41289 KNOTS HEEL AHGLE= 13.76510 DEGREES 


LEEWAY ANGLE= 1.54980 DEGREES RUDDER ANGLE= 1. S9892DEGREES 


7 





SSS SS SSE SSS SESS SSS VSS SEF E Se Se FSO eashneoe Geese @ OSes = _—- 


KYT= 5.00000 KNOTS DGABBA= 85.00000 DEGREES 


KVa= 7.975548 KNOTS DBETTA= 26.31985 DEGREES 


SUCCESSFUL CONVEBGENCE 


ITzZR 


BR 
BOAT SPEED= 3.61350 ‘KNOTS HEEL ASGLZ= 13.57446 DEGREBS 
LEEWAY ANGLE= 2.49007 DEGREES RUDDER ANGLE= 1. 52272DEGREES 
KVT= 5.00000 KsOTS- DGAMMA=  35.00000 DEGREES 
KVA= 7.48743 KNOTS DEBTTA= 22.64893 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


i 
BOAT SPEED= 2.77712 KNOTS HEEL ANGLE= 12.01359 DEGREES 
LEEWAY ANGLE= 3.95441 DEGREES RUDDER AHGLE= 1. 51425DEGREES 
KVT= 7.50000 KNOTS DGABHA= 180.00000 DEGREES 
KVA= 3.76861 KNOTS DEETTA= 179.99980 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


4 
BOAT SPEED= 3.73182 KNOTS HEEL ANGLE= 0.00292 DEGREES 
LEEWAY ANGLE= -0.00395 DEGREES RUDDER ASGLE= 0.01209 DEGREES 
KYT= 7.50000 KNOTS DGABHA=  160.00000 DEGREES 
KVA= 3.97397 KNOTS DEETTA= 139.79790 DEGREES 


SUCCESSFUL CONVERGENCE 


Eee 





ITER = 2 


BOAT SPEED= 4.01511 KNOTS HEEL ANSGLE= 0.05447 
LEEWAY ANGLE= 0.019048 DEGREES RUDDER AGL E= 0.01699DEGREES 
KVT= 7.50000 KNOTS — DG ABB A= 140.00000 DEGREES 


KVA= 4.91665 KNOTS DBETTA= 101.32590 DEGREES 
SUCCESSPUL CONVERGEBRCE 


ITER = a 


BOAT SPEED= 4.77972 KNOTS HEEL ANGLE= 0.86860 

LEEWAY ASGLE= 0.13080 DEGREES RUDDER ANGLE= 0. 0OS872DEGREES 
KVT= 7.50000 KNOTS DG AMS A= 120.00000 DEGREES 

KVA= 6.75196 KNOTS DBETTA= 74.148812 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = 8 


BOAT SPEED= 5.59371 KBOTS HEEL ANGLE= 8.78607 DEGREES 
LEEWAY ASGLE= 0.28889 DEGREES RUDDER ANGLE= 0.56792DEGREES 
KV?= 7.50000 KNOTS DGANBA= 100.00000 DEGREES 

KVA= 8.771484 KNOTS DBETTA= 57.35799 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = 4 


BOAT SPEED= 6.03347 KNOTS HEEL ANGLE= 11.953908 


LEEWAY ANGLE= 0.35957 DEGREES RUDDER AHGLE= 1. 76 328 DEGREES 


ILS) 





KYT= 7.50000 KNOTS DGANMNA= 80.00000 DEGREES 


KYA= 10.36289 KNOTS DPETTA= §5.45833 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


3 
BOAT SPEED= 5.96515 KBOTS HEEL ANGLE= 19.44243 DEGREES 
LEEWAY AWNGLE= 0.50853 DEGREES RUDDER ANGLE= 3. 10435DEGREES 
KYT= 7.50000 KNOTS DGABMA= 60.00000 D2GREES 
KVA= 11.17890 KNOTS DBETTA= 35.52245 DEGREES 


SUCCESSFUL CONVERGENCE 


IT2R 


2 
BOAT SPEED= 5.34839 KNOTS HEEL AWGLE= 24.68011 DEGREES 
LEEWAY ANGLE= 1.17645 DEGREES RUDDER ASGLE= 3.94630DZGREES 
KYT= 7.50000 KNOTS DGASSA= 485.0000C DEGREES 
KVA= 141.17072 KNOTS DPETTA= 28.34309 DEGREES 


SUCCESSPUL CONVERGENCE 


ITlR 


4 
BOAT SPEED= 9.52826 KNOTS HEEL ANGLE= 25.90138 DEGREPS 
LEEWAY ANGLE= 2-448491 DEGREES RUDDER ANGLE= 3.91688 DEGREES 
KYT= 7.50000 KH#OTS DGAMMA= 35.00000 DEGREES 
KVYA= 10.53564 KNOTS DEETTA= 28.09880 DEGREES 


SUCCESSFUL CONVERGENCE 


160 





KvVT= 7.50000 KNOTS DGAMMA= 80.00000 DEGREES 


KVA= 10.36289 KNOTS DEETTA= 45.45833 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


3 
BOAT SPEED= 5.96515 KNOTS HEEL ANGLE= 19.44243 DEGREES 
LEEWAY ANGLE= 0.50853 DEGREES. RUDDER ANGLE= 3. 10935 DEGREES 
KVT= 7.50000 KNOTS DGABBA= 60.00000 DEGREES 
KVA= 11.17890 KNOTS DBETTA= 35.52245 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


4 
BOAT SPEED= 5.34839 KNOTS HEEL ANGLE= 24.68011 DEGREES 
LEEWAY ANGLE= 1.17685 DEGREES RUDDER AEGLE= 3. 94630DZGREES 
KYT= 7.50000 KNOTS DGABBA= &5.0000C DEGREES 

KVA= 11.17072 KNOTS DEETTA= 28.34309 DEGREES 


SUCCESSPUL CONVERGENCE 


ITZR 


4 
BOAT SPEED= 8.52826 KNOTS HEEL ANGLE= 25.90138 DEGREES 
LEEWAY ANGLE= 2-49491 DEGREES RUDDER ANGLE= 3.91688 DEGREES 
KYT= 7.50000 KsOTS DGAMAA= 35.00000 DEGREES 
KVA= 10.53569 KNOTS DEETTA= 28.09880 DEGREES 


SUCCESSPUL CONVERGENCE 
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ITER 


8 
BOAT SPEED= 3.47308 KNOTS HEEL ANGLE= 23.64365 DEGREES 
LEEWAY ANGLE= &.47674 DEGREES RUDDER AHGLE= 4. OV634 DEGREES 
KVT= 10.00000 KNOTS | DG AMA A= 180.00000 DEGREES 
KV¥A= 5.16953 KNOTS DBETTA= 179.99990 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


& 
BOAT SPEED= 4.83037 KNOTS HEEL ANGLE= 0.00592 DEGREES 
LEEWAY AWGLE= -0.00401 DEGREES RUDDER ANGLE= 0. 01243 DEGREES 
KVT= 10.00000 KNOTS DGAMMA= 160.00000 DEGREES 
KVA= 5.38727 KNOTS DBETTA= 141.84290 DEGREES 


SUCCESSFUL CONVERGENCE 


ITZR 


3 
BOAT SPEED= 5.10578 KNOTS HEEL ANGLE= 0.08175 DEGREES 
LEEWAY ANGLE= 0.01586 DEGREES RUDDER ANGLE= 0. 0277S5DEGREES 
KVT= 10,.00000 KNOTS DG AAMA= 180,.00000 DEGREES 
KVA= 6.67725 KNOTS DBETTA= 105.70840 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER 


3 
BOAT SPEED= 5.85245 KHOTS HEEL ANGLE= 1.40812 DEGREES 
LESWAY ANGLE= 0.10124 DEGREES RUDDER ANGLE= 0. 17259 DEGREES 
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KVYT= 10.00000 KNOTS DGABBA= 120.00000 DEGREES 
KVA= 8.79842 KNOTS DBETTA= 79.97887 DEGREES 
SUCCESSFUL CONVERGENCE 


Itzk = 8 


BOAT SPEED= 6.52654 ‘KNOTS HEEL ANGLE= 6.55590 DEGREES 
LEEWAY ANGLE= 0.10446 DEGREES RUDDER AWGLE= 1. 10168DEGREES 
KVT= 10.00000 KNOTS. DGAMMA=  100.00000 DEGREES 

KVA= 10.96851 KNOTS DBETTA= 63.87666 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 2 § 


BOAT SPEED= 6.56510 KNOTS HEEL ANGLE= 16.25755 DEGREES 
LEEWAY ANGLE= -0.16120 DEGREES RUDDER AXGLE= 3.55508 DEGREES 
KVT= 10.00000 kNOTS DGAMAA= 80.00000 DEGREES 

KVA= 12.85525 KNOTS DBETTA= 50.00259 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = 3 


BOAT SPEED= 6.52192 KNOTS HEEL ANGLE= 27.24298 DEGREES 
LEEWAY ANGLE= -0.33663 DEGREES RUDDER ANGLE= 5. 98752DEGREES 
KVT= 10.00000 KNOTS DGAMBMA= 60.00000 DEGREES 

KVA= 13.80896 KNOTS DBETTA= 38.84161 DEGREES 


SUCCESSFUL CONVERG EWCE 


G2 





ITER = 5 


BOAT SPEED= S.754814 KNOTS HEEL ASGLE= 35.06798 DEGREES 
LEEWAY ANGLE= 0.37846 DEGREES RUDDER ANGLE= 7. 6084S5DEGREES 
K VT= 10.00000 KNOTS . DGASSA= 45.00000 DEGREES 


KVA= 13.68682 KNOTS DBETTA= 31.10666 DEGREES 
SUCCESSPUL CONVERGESCE 


ITER = & 


BOAT SPEED= 4.64757: KSOTS HEEL ANGLE= 37.00015 DEGREES 
LEEWAY ANGLE= 2.45695 DEGREES RUDDES ANGLE= 8, 19907 DEGREES 
K VT= 10.00000 KNOTS DGABHA= 35.00000 DEGREES 

KVA= 12.89910 KNOTS DBETTA= 27.31577 DEGREES 


SUCCESSPUL CONVERGENCE? 
ITER = 7 
BOAT SPEED= 2.91379 KNOTS HEEL ANGLE= 32.23439 DEGREES 


LEZWAY ANGLE= 7.24246 DEGREES RUDDER ANGLE= 13. 07657 DEGREES 
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ITMAX = 10 
IPRINT = 0 
y « rf 
EPSi1 = 1. OD- 10 
EPS2 =, 1. 0B-02 


XOLD(1)...XOLD( 4) 


1,0000002 01 0.0000002-01  0.0000002-01 0.0000008-01 
SOLOUTIOHS ARE COMPUTED FOR THE FOLLOWING TRUE WIND SPEEDS (KNOTS) 
§.00000 7.50000 10.00000 


SOLUTIONS ARB COMPUTED FOR THE FOLLOWING TRUE WIND ANGLES FOR EACH WIND SPEED 


180.00000 160.00000 140.00000 120.00000 100.00000 80.00000 60.00000 45.00000 35.00000 


» 
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LWL= 30.00000 

DISPLACEMEST TO LENGTH*  225.00000 
PRIS SATIC= 0.53000 

L/a= 14. 76796 


LCB= 53.00000 
KIHVEMATIC VISCOSITY= 0.1279100E-038 


DENSITY OF WATER= 1.99000 

MEAN KEEL CHORD= 6.90000 

MEAN RUDDER CHORD= 1.87537 ‘ 
ROUGHNESS AND CORRELATION FACTOR= 0.00000 


PAIRBODY WETTED SURPACE= 215.25700 

KEEL WETTED SURPACE= 63. 01086 

RUDDER WETTED SURFACE= 19.99725 

TOTAL WETTED SURPACE= 298.26480 

BYL /L&L= 0.30276 

SAIL AREA=  900.00000 

HEELING MCMENT ARA= 23.99541 

PITCHING MOMENT ARA= 26. 35417 

DISTANCE FROM FP TO RUDDER STOCK= 28. 49998 
CLOSE HAULED LIPT COEP= 1.25000 

RIG PORA DRAG COEF= 0. 90000 

CROSS SECTIONAL AREA OF RIG= 26.97595 

SAIL PRICTIONAL COEF= 0.09000 

INDUCED DRAG COZP / LIFT COEP**2= 0. 13030 


LATERAL AREA OF HULL= 121. 29310 


TRANSVERSE AREA OF HULL= 37.06689 
DISTANCE BETWEEN CE SAIL AND PP= 13.93835 
OPP BIND DRAG COEF= 0.00000 
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OFF WIND LIFT COEP= 0.00000 

DISTANCE BETWEEN CE HULL AND PP= 13.93835 

DENSITY OF AIR= 0.00238 

HORIZONTAL DISTANCE PROM MAXISUM DRAFT OF PROFILE MEASUREMENTS 
0.00000 3.00000 6.00000 9.00000 16.50000 

PROPILE MEASUREMENTS 


6.38000 “6.25240 1.92985 1.42200 0.00000 


DRAFT OF CANOE BODY / TOTAL DRAFT= 0.31841 
‘ 
HYDRODYNAMIC LIFT COEP=— 0.11209 
LANDA X LIFT / DRAG= 2.08978 
DISTANCE FROM LWL TO VERTICAL HYDRO CE= 2.12954 


*** RUDDER SECTION *** 


DISTANCE PROM LWL TO ROOT CHORD= 0.40628 
ROOT CHORD= 2.25000 

TIP CHCRD= 1.50075 

RUDDER SPAN= 5.07766 

SWEEP ANGLE= 0.00000 

ASPECT RATIO= 5.41509 

MEAN CHORD= 1.87537 

DCL / DALPHA= 0.07057 

CROSS PLOW DRAG COEFP= 0.56024 
DCS / DCL= 0.21654 

FORA DRAG COEFP= 0.00650 
RUDDER AREA= 9.52251 


*** RIGHTING MOMENT SECTION *** 
METACENTRIC HEIGHT / LWL= 0.11534 


FROUDE NUMBER RIGHTING MOMENT COEF= 0.60895 152-04 
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BON LIHEAR HYDROSTATIC COEF= 0.3767873 2-03 


DISTANCE BELOW LWL TO CG= 0.75000 
**e SAIL SECTION #8 


I= 47. 17590 d=: 18.72513 

P= 82. 88670 " Be 14.29557 

AB POR TRI= 3.00000 AR SAIN= 3.00000 
AB SAIL PLAW= 2.47280 

CE SAIL PLAN APT MAST= 0.60041 

FORWARD OVERHANG= 2.38720 

F REE BOARD = 3.71000 


ILS) 7 
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KYT= 5.00000 KNOTS DGABASA= 180.00000 DEGREES 
KVA= 2.36893 KNOTS DBETTA= 179.99980 DEGREES 
SUCCESSFUL CONVERGENCE 


ITER * S 


BOAT SPEED= 2.63557 KNOTS HEEL ABGLE= 0.00057 DEGREES 
LEEWAY ANGLE= -0.00116 DEGREES RUDDER ANGLE= 0. 00433DEGREES 
KYT= 5.00000 KNOTS DGARMOA= 160.00000 DEGREES 

KVA= 2.49139 KNOTS DBETTA= 136.65290 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER = 2 


BOAT SPEED= 2.88579 KNOTS REEL ANGLE= 0.01554 DEGREES 


LEEWAY ANGLE= 0.02C97 DEGREES RUDDER ANGLE= 0. 00262DEGREES 


WO CONVERGENCE 
KYT= 5.00000 KNOTS DGAMMA= 140.00000 DEGREES 


KVA= 3.21880 KxOTS DBETTA= 93.15175 DEGREZS 
SUCCESSFUL CONVERGESCE 


ITER = 11 


BOAT SPEED= 3.65326 KNOTS HEEL ANGLE= 0.24055 DEGREES 
LEEWAY ANGLE= 0.13603 DEGREES RUDDER ANGLE= 0, 02736 DEGREES 
KVt= 5.00000 KNOTS DGANSA=  120.00000 DEGREES 

KVA= 4.83768 KNOTS DBETTA= 63.51920 DEGREES 


SUCCESSFUL CONVERGENCE 
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ITER 


8 
BOAT SPEED= 4.65673 KNOTS HEEL ANGLE= 1.80063 DEGREES 
LEEWAY ANGLE= 0.54518 DEGREES RUDDER ANGLE= 0. 06240 DEGREES 
KYT= 5.00000 KsOTS DGAMBA= 100.00000 DEGREES 


KVA= 6.53116 KNOTS DBETTA= 48.93198 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


& 
BOAT SPEED= 5.15885 KNOTS HEEL ANGLE= 4.21832 DEGREES 
LEEWAY ANGLE= 0.98065 DEGREES RUDDER ANGLE= 0. 20861DEGREES 
KVT= 5.00000 KNOTS DGAAMA= 80.00000 DEGREES 
KVA= 7.84539 KWOTS DBETTA= 38.87585 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER 


3 
BOAT SPEED= 5.23940 KNOTS HEEL ANGLE= 6.78553 DEGREES 
LEEWAY ANGLE= 1.49549 DEGREES RUDDER ANGLE= 0. 37304DEGREES 
KVT= 5.00000 KNOTS DGABAA= 60.00000 DFPGREES 
KVA= 8.57841 KNOTS DBETTA= 30.316%12 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


3 
BOAT SPEED= %.90367 KNOTS HEEL ANGLE= 8.62737 DEGREES 
LEEWAY ANGLE= 2.718603 DEGREES RUDDER ANGLE= 0. 45298DEGREES 
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KYT= 5.00000 kKwOoTS DGARBA= 85.00000 DEGREES 
KVA= 8.39471 KNOTS DEETTA= 28.90797 DEGREES 
SOCCESSPUL CONVERGENCE 


ITER = 4 


BOAT SPEED= 8.07831 KxOTS HEEL ANGLE= 8.89245 DEGREES 
LEEWAY ABNGLE= 3.16383 DEGREES ; RUDDER ANGLE= 0. 93058 DEGREES 
KVT= 5.00000 KwWOTS DGAMMA= 35.00000 DEGREES 

KVA= 7.81175 KgOTS DBETTA= 21.53836 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = nh 


BOAT SPEED= 3.17022 KNOTS HEEL ANGLE= 7.84635 DEGREES 
LEEWAY ANGLE= 4.64773 DEGREES RUDDER ASGLE= 0. 38519 DEGREES 
KYT= 7.50000 KNOTS DGAMMA=  180.00000 DEGREES 

KYA= 3.59100 KNOTS DBETTA= 179.99980 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = 3 


BOAT SPEED= 3.91276 KHOTS HEEL ANGLE= 0.00121 DEGREES 
LEEWAY ANGLE= 0.00031 DEGREES RUDDER ANGLE= 0. 00054 DEGREES 
KYT= 7.50000 KKOTS DGABMA=  160.00000 DEGREES 

KYA= 3.81778 KHOTS DBETTA= 137.78660 DEGREES 


SUCCPSSPUL CONVERGENCE 


IL 718, 





ITER 


2: 


BOAT SPEED= 4.21828 KNOTS HEEL ABGLE= 0.03173 DEGREES 
LEEWAY ANGLE= = 0.02009 DEGREES RUDDER ANGLE= 0. 00629DEGREES 
KVT= 7.50000 KwoTS: DGAMHA=  190.00000 DEGREES 


KVA= 4.86734 KNOTS DBETTA= 97.92107 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER 


3 
BOAT SPEED 5.07021 KNOTS HEEL ANGLE= 0.89961 DEGREES 
LEEWAY ANGLE= 0.14829 DEGREES RUDDER ANGLE= 0. O1619DEGREES 
KVT= 7.50000 KNOTS DGASHMA= 120.00000 DEGREES 
KVA= 6.89121 KNOTS DBETTA= 71.699138 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER 


3 
BOAT SPEED= 5.90355 KmsoTS HEEL ANGLE= 2.91928 DEGREES 
LEEWAY ABGLE= 0.51083 DEGREES RUDDER ANGLE= 0. 18288DEGREES 
KYT= 7.50000 KNOTS DGAHBSA= 100.00000 DEGREES 
KVA= 9.05576 KNOTS DBETTA= 54.64867 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER 


3 
BOAT SPEED 6.53885 KNOTS HEEL ANGLE= 7.42339 DEGREES 
LEEWAY ANGLE= 0.95508 DEGREES RUDDER AUGLE= 0. 56130DEGREES 


SSX Q eae ee Se Seo Sl 8S SBS] 2246 B28 SF S228 SF SEBS SS 2S SOS SCF O82 BSH 


1 ag 





K ¥T= 7.50000 KsOTS DGABMA= 80.00000 DEGREES 
KYVA= 10.80911 Ks8OTS DBETTA= 43.10342 DEGREES 
SUCCESSFUL CONVERGENCE 
ITER = 3 


BOAT SPPED= 6.58876 KNOTS HEEL ANGLE= 12.24898 DEGREES 


LEPUAY ANGLE= 1.51378 DEGREES : RUDDER ANGLE= 0. 99877 DEGREES 
KVT= 7.50000 KNOTS’ DGABNA= 60.00000 DEGREES 
KVA= 11.78616 KNOTS DBETTA= 33.57051 DEGREES 


SUCCESSPOL CONVERGENCE 


ITER = & 


BOAT SPEED= 6.03699 KHOTS BEEL ANGLE= 15.686428 DEGREES 
LEPWAY ANGLE= 2-41843 DEGREES RUDDER ABGLE= 1, 26 121DEGREES 
KVT= 7.50000 KNOTS DGABMA= &85.00000 DEGREES 

KVA= 11.81018 Kxrots DBETTA= 26.68242 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER = 3 


BOAT SPEED= 5.24840 KNOTS HEEL ABGLE= 16.53915 DEGREES 
LEEWAY ANGL2= 3.57327 DEGREES RUDDER ANGLE= 1, 20193DEGREPS 
KYT= 7.50000 KNOTS DGARSA= 35.00000 DEGREES 

KVA= 11.24568 KNOTS DBETTA= 22.89055 DEGREES 


SUCCESSFUL COBVERGENCE 


NZ 





ITER 


a 


SUCCESSFUL CONVERGENCE 


ITER 


SUCCESSPUL CONVERGENCE 


ITER 


SUCCESSFUL CONVERGENCE 


ITER 


BOAT SPEED= * 8.24629 KNOTS HEEL ANGLE= 15.29806 DEGREES 
LEEWAY ANGLE= 5.22368 DEGREES RUDDER ANGLE= 1. O9039DEGREES 
KvT= 10.00000 KNOTS , DGAMNM A= 180.00000 DEGREES 
KVA= 9.96443 KNOTS DBETTA= 179.99980 DEGREES 

3 
BOAT SPEED §.03558° KNOTS HEEL ANGLE= 0.00228 DEGREES 
LEEWAY AHGLE= -0.0004%6 DEGREES RUDDER ANGLE= 0.00276DEGREES 
KvT= 10.00000 kKSOTS DGAMMA= 160.00000 DEGREES 
KVA= 5.32662 KyOTS DBETTA= 180.05190 DEGREES 

3 
BOAT SPPED= 5.31337 KNOTS HEEL ANGLE= 0.04492 DPGREES 
LEEWAY ANGLE= 0.01862 DEGREES RUDDER ANGLP= 0.00731DEGREES 
KVT= 10.00000 KwOoTS DGAREA= 180.00000 DEGREES 
KVA= 6.61337 KNOTS DBETTA= 103.60230 DPGREES 

3 
BOAT SPPED= 6.10511 KNOTS BEEL ANGLE= 0.78191 DEGREES 
LEEWAY ANGLE= 0.14843 DEGREES RUDDER ANGLE= 0.04948 DEGREES 


L73 





KV¥T= 10.00000 KxoTS DGAMAA= 120.00000 DEGREES 


KVA= 8.88892 KNOTS DBETTA= 77.08757 DEGREES 


SUCCESSPUL CONVERGENCE 


ITER 


9 
BOAT SPEED= 6.98967 KNOTS HEEL ANGLE= 4.09941 DEGREES 
LEEWAY ANGLE= 0.45382 DEGREES RUDDER ANGLE= 0. 3520S5DEGREFS 
KYT= 10.00000 KsoTS DGANMA=  100.00000 DEGREES 
KVA= 11.29239 KNOTS DBETTA= 60.70335 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


8 
BOAT SPEED= 7.26006 KNOTS HEEL ANGLE= 10.27948 DEGREES 
LEEWAY ANGLE= 0.89417 DEGREES RUDDER ANGLE= 1. 12386 DEGREES 
KYT= 10.00000 xKsoTs DGABMA= 80.00000 DEGREES 
KVA= 13.33439 KNOTS DBETTA= 47.60773 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER 


3 
BOAT SPEED= 7.25288 KNOTS HEEL ASGLE= 17.39751 DEGREES 
LEEWAY ANGLE= 1.48165 DEGREES RUDDER ANGLE= 2. 01724 DEGREES 
KVT= 10.00000 KNOTS DGABMA= 60.00000 DEGREES 
KVA= 14.65839 KNOTS DBETTA= 36.21914% DEGREES 


SUCCESSFUL CONVERGENCE 
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ITER = 5 


BOAT SPEED= 6.82657 KNOTS HEEL ANGLE= 23.94096 DEGREES 
LEEWAY ANGLE= 2-54576 DEGREES RUDDER ANGLE= 2. SO608 DEGREES 
K VT= 10.00000 KROTS DGAMBA= 45.00000 DEGREES 

KVA= 14.70852 KNOTS DBETTA= 28.74271 DEGREES 


SUCCESSFUL CONVERGENCE 


IT 28 = 4 


BOAT SPZED= 5.82167 KNOTS HEEL ANGLE= 24.95070 

LEEWAY ABGLE= 4.21530 DEGREES RUDDER ANGLE= 2.51386DEGREES 
KVT= 10.00000 KNOTS DGAMBA= 35.00000 DEGREES 

KVA= 18.17157 KNOTS DEETTA= 23.879566 DEGREES 


SUCCESSFUL CONVERGENCE 


ITER = & 


BOAT SPEED= 4.76666 KNOTS HEEL ANGLE= 23.88383 


LEEWAY ANGLE= 6.43894 DEGREES RUDDER ANGLE= 2. 34001 DEGREES 
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